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FOREWORD 


This  report  was  prepared  by  Aeronautical  Research 
Associates  of  Princeton,  Inc.,  Princeton,  New  Jersey  under 
Air  Force  Contract  AF33( 615 )-2878,BPSN: 5 (61136702-62405334). 

The  contract  was  Initiated  under  Project  No.  1367,  "Structural 
Design  Criteria,  Task  No.  136702,  "Aerospace  Vehicle  Structural 
Loads  Criteria". 

The  work  was  administered  under  the  direction  of  the  Air 
Force  Flight  Dynamics  Laboratory,  Research  and  Technology 
Division,  Air  Force  Systems  Command,  Wright- Patter son  Air 
Force  Base,  Ohio,  Mr.  Paul  L.  Hasty  (FDTR),  Project  Engineer. 

The  work  reported  in  this  study  was  conducted  by  Aeronau¬ 
tical  Research  Associates  of  Princeton,  Inc.  with  Dr.  John 
C.  Houbolt  as  principal  Investigator,  and  covers  the  period 
April  1966  to  April  1967 . 

This  technical  report  has  been  reviewed  and  is  approved. 
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ABSTRACT 


Further  developments  of  gust  design  procedures  based  on 
power  spectral  techniques  are  presented-  A  number  ol  topics 
are  considered,,  and  presentation  is  in  the  nature  of  a  series 
of  interrelated  small  reports  under  ohe  cover.  Generalized 
load  exceedance  curves  are  considered  theoretically,  and 
certain  significant  properties  of  these  curves  are  established. 
It  is  shown  that  response  and  design  considerations  may  be 
expressed  in  terms  of  three  basic  parameters :  P  the  propor¬ 
tion  of  time  in  turbulence,  a_  the  gust  severity,  and  a  shape 

v 

parameter  which  defines  the  generalized  exceedance  curves.  Over 
a  dogen  different  families  cf  theoretical  exceedance  curves 
are  generated.  The  basic  response  and  environmerfc al  parameters 
that  are  of  concern  in  design  are  discussed,  and  the  use. of 
composite  values  of  these  parameters  as  might  be  involved  in 
mission  considerations  is  shown.  The  composite  gust  Intensity 
value  e?c  ,  and  the  related  scale  value,  L  ,  still  represent 

unsettled  questions.  A  reexamination  of  seme  previous  airline 
operational  data  in  terms  of  generalized  exceedance  curves  is 
included . 

% 

Recommendations  on  four  specific  design  procedures  are 
given.  One  of  the  design  procedures,  based  primarily  on  the 
response  parameters  A  and  ,  incorporates  the  results  of 

computational  studies  that  were  performed  on  certain  existing 
aircraft  as  a  means  for  establishing  design  boundaries.  -Rie 
problem  of  determining  the  probability  of  exceeding  given  load 
levels  in  flights  of  specified  duration  is  investigated.  Areas 
and  parameters  which  are  considered  to  have  weakness  or  uncer¬ 
tainty  are  indicated,  and  recommendations  are  accordingly  made 
for  appropriate  future  research  effort.  Consideration  of  the 
effect  of  filtering  on  deduced  scale  value  is  included  in  an 
Appendix. 


This  abstract  is  subject  to  specific  export  controls,  and 
each  transmittal  to  foreign  governments  or  foreign  nationals 
may  be  made  only  with  prior  approval  of  the  Air  Force  Flight 
Dynamics  Laboratory  (FDTR),  #righfc-Patterson  Air  Fores  Base, 
Ohio  45433 • 
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IHTEQBOCnOg 


Oncer  contract  with  the  Air  Force  Flight  Dynamics  Labora¬ 
tory,  Aeronautical  Research  Associates  of  Princeton  has  been 
developing  procs ceres  for  designing  aircraft  for  gast  encounter 
based  on  power  spectral  techniques.  Reference  1  reports  the 
prelistinary  findings  of  these  st Tidies.  In  particular,  four 
possible  design  procedures  covering  various  degrees  of  sophistica¬ 
tion  were  brought  out.  Sirs  report  presents  the  results  of  a 
con  tinning  phase  of  effort  aimed  at  evolving  further  the  proced¬ 
ures,  as  well  as  investigating  other  aspects  of  the  spectral 
approach  to  gust  response  and  design  problems.  The  presenta¬ 
tion  given  herein  is  in  the  nature  of  a  series  of  small  reports, 
covering  a  variety  of  interrelated  hspics  which  pertain  to  the 
subject  of  gust  loads  analysis  of  aircraft. 

As  W2S  brought  oat  in  Reference  1,  related  work  has  simul¬ 
taneously  been  carried  out  by  some  of  the  aircraft  companies, 
sponsored  also  by  the  Air  Force  Flight  Dynamics  Laboratory 
and  under  Ron!  tor  ship  by  Dr.  J.  C.  Eoubolt,  the  principal  investi¬ 
gator  of  the  present  investigation.  These  were  computational 
studies  of  certain  existing  gust  critical  airplanes  that  were  *  ! 

designed  by  tbs  discrete-gust  technique,  and  were  aimed  at  j 

establishing  for  these  aircraft  some  of  the  basic  response 
parameters  that  are  significant  in  a  power  spectral  design  approach. 
The  results.  References  Sk.  9,  and  10,  are  incorporated  herein 
as  one  of  the  means  for  guiding  the  choice  of  spectral  design 
numbers  and  design  boundaries.  It  is  well  to  repeat  here  one 
of  the  underlying  theses  that  is  being  followed,  namely,  that  of 
applying  a  possible  new  procedure  in  retrospect  manner  to  proven 
gust  critical  aircraft  with  the  concept  that  if  they  had  been 
designed  originally  by  the  new  procedure,  safe  designs  would 
have  also  resulted. 


SECTION  2 


THEORY  FOR  IEVBL  CROSSINGS,  Nfx) 


This  section  gives  some  of  the  basic  theory  that  is 
associated  with  level  or  threshold  crossings  for  random  time 
history  functions  «.  Application  is  then  made  to  the  gust  leads 
case,  'with  the  view  toward  developing  further  the  concept  of 
nuniversal  load  exceedance  curves’1  that  is  brought  out  in 
Reference  1.  It  is  shown,  in  quite  general  form,  that  load 
exceedance  considerations  need  only  involve  three  basic 
environmental  factors.  A  new  and  remarkable  moment  property 
of  the  level  crossing  curves  is  .also  established  which  is  of 
governing'  importance  in  the  mathematical  description  of  these 
curves  and  which  should  he  of  marked  significance  in  the 
interpretation  and  use  of  flight  data  on  load  exceedance. 

r 

General  Formulation.-  A  derivation  of  the  general  equa¬ 
tion  for  determining  the  number  of  level  crossings  is  given 
first,  using  an  approach  which  is  slightly  different  than  that 
used  -by Uice  in  Reference.  2,.  Consider  a  random  function  of 
time  as  depicted'  in  the  following  sketch: 


dx 


We  wish  ip  establish  the  statistical  estimate  for  the  average 
number  of  crossings  of  the  level  x  .  For  any  one  crossing, 
the ... amount  -of  time  to  cross  the  interval  dx  is 


In  a  total  time  T  ,  the  total  time  spent  in  the  incremental 
area  dxdx  is,  by  definition 


AT  =  Tp(x,x)dxdx 


where  p  represents  the  joint-probability  distribution  function 
of  the  random  function  x  and  its  time  derivative  x  .  The 
total  number  of  dx  crossings  in  time  T  is  thus 
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dn  =  ^  -  Tjxjp(x,x)dx 


We  now  find  the  contribution  from  all  possible  x‘s  by  inte¬ 
grating  with  respect  to  x  *  thus. 

CO 

n  =  T  J  jxjp(x*x)dx  ' 

—CO 


The  probability  density  function  is  supposed  symmetric*  and* 
therefore*  the  -  contribution  from  negative  x  should  be  the 
same  as  for  positive  x  ;  the  equation  may  tnus  be  written 

CO 

n  =  2Tjf  x  p(x*x)dx  (1) 

0 

This  equation  yields  the:  total  number  of  crossings,  both  up¬ 
ward  and  downward*  of  the  level  x  .  The  number  of  upward 
crossings  (or  the  number  of  downward  crossings  for  negative 
x)  is  given  with  the  factor  2  suppressed.  Thus*  suppressing 
the  2  and  putting  the  equation  on  a  per  unit  time  basis,  we 
obtain  the  following- general  form  for  the  number  of  upward 
crossings  per  unit  time  of  the  level  x 

00 

M(x)  =f  =/xp(xJx)dx  (2) 

0 


This  equation  is  true  regardless  of  the  form  of  p.  Specific 
analytical  egressions  for  N(x)  may  be  found  by  substituting 
in  any  plausible  and  convenient  expression  for  the  joint 
distribution  function  p  .  In  general*  these  joint  distri¬ 
bution  functions  need  only  have  the  following  constraints 
or  properties; 


00  00 


1  -  J'  J  p(x,x)dxdx 

—00  -00 
00  00 

°x2  =//x2P(  x*x)dxdx 
—00  —00 
CO  00 

Cjj*  = x^p(x*x)dxdx 


(3) 


In  addition,  we  limit  ourselves  to  only  those  p  functions 
which  are  symmetrical  in  "both  the  x  and  x  directions. 

A  study  of  these  equations  indicates  that  p{x,x)  must 
be  of  the  form 

p(x^}  k’ a)  m 


where  g  ia  a  nondimensionai  function  which  is  described 
completely  by  means  of  a  “shape”  parameter  a  ,  and  by  the 

nondimensionai  independent  variables  x/a  and  x/cr*  , 

^  ^  « 

which  involve  a  and  a*  ,  the  r.m.s.  values  of  x  and  x  . 

X  X 

The  factor  l/2w  Is  included  simply  as  an  aid  to  be  used  later 
in  defining  the  basic  frequency  parameter  NQ  .  (Note  the 

symmetric  properties  of  p  can  be  explicitly  brought  out  by 

iza 

simplicity  the  square  notation  will  be  omitted.)  Substitution 
of  this  equation  for  p  into  Equations  3,  and  introducing 
the  notation 


^  but  for 


writing  g  in  the  functional  form 


x 


(5a) 


n  = 


X 


(5b) 


gives  the  following  properties  of  g 


00  00 


(6) 


(7) 
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1  =  hr  f  /n2eUi  n>  <*) 


-CO  -00 


(8) 
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p(*) 


CO 


1_  1_ 
2ir  u 


f  &(Z,  cl) 


dT) 


(9) 


Application  to  Gust  Encounter  -  The  Three  Basic  Environ¬ 
mental  Factors  *  P  ,  and  a  .-  in  application  to  the 

airplane  gust  encounter  problem*  Equation  2  applies*  of  course* 
only  during  turbulence  encounter.  To  make  the  equation  yield 
results  which  give  the  average  level  crossing  rate  based  on 
entire  flight  time*  we  multiply  by  T^  *  the  total  time  spent 

in  turbulence*  and  divide  by  T  *  the  total  flight  time:  the 
result  is 

CO 

N(x)  =  P  jf  xp(x,x)dx  (10) 

0 

where  P  =  T^/T  denotes  the  proportion  of  time  spent  in  tur¬ 
bulence.  Substitution  of  Equation  4  into  this  equation  yields 


00 


This  equation  is  of  fundamental  signif  cance.  Apart  from  * 

which  is  related  to  the  structural  response  characteristics  of 
the  aircraft*  it  shows  that  the  level  crossing  history  of  an 
aircraft  in  gust  is  dependent  on  only  three  basic  parameters: 

P  *  the  proportion  of  time  in  turbulence;  a x  *  the  severity 

parameter;  and  a  the  shape  parameter  of  the  curve  f  ' )  . 

It  is  noted  that  the  equation  is  derived  through  general  con¬ 
sideration  only.  Concepts  of  an  atmospheric  model  consisting 
of  a  series  of  patches  (see  Reference  3)*  or  Involving  the 
assumption  of  Gaussian  distribution,  need  not  be  introduced. 
Instead  the  concept  is  advanced  that  the  gust  encounter  ex¬ 
perience  of  an  aircraft  can  be  described  completely  in  terms 
of  an  overall  distribution  curve  of  the  type  given  by  Equation 
4.  or*  alternatively*  by  Equation  11.  These  two  equations 
form  the  basis  for  developments  contained  in  this  report. 
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Sees  general  observations  pertaining  to  the  ^slieahility 
of  Equation  11  sdgbt  sects  at  this  point.,  Firsts  x 
represents  any  variable  of  concern,  either  input  or  curtpst 
response.  If  f  is  a  response  variable.,  then  its  r .a. s . 
value  is  related  to  the  r.m.s.  gust  value  through  the  weil- 
•down  expression 


The  "severity11  c,r  is  thus  related  directly  to  the  composite 

severity  oc  of  the  gusts  through  the  parameter  A  .  "The 

practical  evaluation  of  the  structural  response  quantities  A 
and  Hq  sy  spectral  techniques  is  considered  in  “Reference  1 

and  3. 

Secondly,  we  note,  without  detailed  elaboration,  that 
ecirpleia  stationarity  is  not  a  necessary  prerequisite  for  the 
randan  variable  x  .  Stationarity  of  is  Implied  but  the 

local  r.m.s.  value  (as  defined  on  a  patch  sense  basis)  need 
not  be.  Varying  HQ  values  are  easily  taken  into  account  by 

a  superposition  technique. 

Thirdly,  it  may  be  remarked  that  the  observations  sug¬ 
gested  by  Equation  11  are  also  supported  by  practical  consi¬ 
derations.  Thus,  if  one  reflects  and  tries  to  single  out  the 
main  factors  which  govern  gust  loads  experience,  it  seems 
plausible  that  the  logic  would  narrow  down  to  the  following 
three  questions: 

1)  What  is  the  time  spent  in  turbulence? 

2)  What  is  the  severity? 

3)  What  is  the  general  shape  of  any  chosen 
load-describing  curve? 

These  three  questions  are  reflected  by  Equation  .11. 

In  the  next  section  it  will  be  shown  that  the  first  mo¬ 
ment  of  the  area  under  the  N(x)  curve  about  the  origin  (con¬ 
sidering  the  right-hand  plane  only)  is  of  significant  practi¬ 
cal  importance .  Thus  with  Equations  4  and  10  we  write  the 
following  for  later  use 


(12) 


•7 


General  Besults  and  Basic  Kogent  Property  For  a  Practical 
Subclass If  we  restrict  oar  attention  to  a  certain  practical 
class  of  randan.  functions  where*  at  least*  the  function  and  the 
first  derivative  are  continuous*  then  further  properties  of 
note  say  he  deduced  in  general .  Specifically,  suppose  the 
fractions  are  described  by  joint -distribution  functions  wherein 
the  independent  variables  always  appear  in  the  combination 


=  r 


o 

x 


5Te  introduce  the  change  of  variable 


e  =  - 

5  *x 


=  — —  =  r  cos  6 


n  —  — — -  =  r  sin  6 
ax 


Then  for  a  single  integration  in  the  ij  direction 


T]  =  A2  - 


iti  = 


rdr 


r~2.  7i 

V  T*  -  t 


and  for  double  integration*  we  have*  through  the  Jacobian 


d£dt]  =  rdrd9 

With  these  relations*  Equations  6*  7>  11*  and  12  become 

00 

r\ 

.1  =  J  rg(r,  a)  dr 


o 


From  the  first  and  third  of  these  equations,  we  see  that 

or  f  (x,u)  of  Equation  11  must  have  the  property 

- 

=  t  (0,a)  =  1  (13) 


while  the  second  and  fourth  indicate  that 


00 


Equation  14  is  of  fundamental  Interest  and  significance. 
Geometrically  it  indicates  that  the  first  moment  of  the  area 
under  the  right-hand  side  of  the  normalized  exceedance  curve 
2 

must  equal  a  ,  a  fact  evidently  not  known  before.  This 

fact  and  Equation  13  place  important  restraints  on  the  form 
that  the  exceedance  curve  can  take,  and  are  most  useful  in 
establishing  analytical  expressions  for  the  exceedance  curve, 
as  will  he  seen.  These  two  properties  should  also  he  very 
helpful  in  the  interpretation  of  experimental  data;  in  fact, 
if  in  Mi e  deduction  of  exceedance  curves  from  flight  data  it 
is  found  that  these  properties  are  not  met,  particularly  the 
moment  property,  then  it  may  be  surmised  that  there  is  some¬ 
thing  wrong  with,  the  data  evaluation  or  interpretation. 

It  is  noted  that  Equation  14  is  exact  if  is  invariant 
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SECTION  3 

SPECIFIC  LOAD  EXC2EDMCE  CURVES 


Derivation  Procedures.-  At  least  three  different 
procedures  may  he  used  for  "deriving  specific  load  exceedance 
curves  as  described  in  the  following. 

Direct  means: 

In  this  procedure,  a  form  of  p(x,x)  satisfying  Equa¬ 
tions  6,  7a  and  8  is  assumed.  Substitution  into  Equation 
10  leads  then  directly  to  the  result  for  N{x)  .  As  an  exam¬ 
ple^  a  cade  often  quoted  and  used  is  "based  on  the  assumption, 
that  x;  '.and  x  obey  a  .joint -normal  or  Gaussian  distribution. 
In  this  instance*  p  is  given  by 


p(x*x)  = 


1_  1 

27r  Vx 


Substitution  of  this  equation  into  Equation  10  and  the  last 
of  Equations  3  yields 


C-U 

N(x).=  PN_e  x 


P(x)  =  ~i_  A-  e  x 


which  are  noted  to  be  a  form  of  Rice’s  results. 


It  is  significant  to  note  that  a  Gaussian  assumption  is 
not  essential  in  the  treatment,  of  gust  loads,  as  has  been  the 
general  belief  or  the  underlying  assumption  in  previous  studies. 
Any  plausible  form  of  p(x?x)  may  be  used.  As  a  specific 
example if 

a 


*(*’*)  =  2 707 


r2T^r. 


^)] 


-  *  . 

5-  .%£**’■ 


3-  ?  * '  ” 


then  the  following  equation  results  for  2J(x) 


I&L  = 


By-  suiifetidh  of  Russian  ease's* 

’  Although  not' neesssaryj  the  technique  described  in 
Reference  1  for  the  composite  gust  model  consisting  of  the 
encounter  of  a-  nu mber  of  Gaussian  patches  of  turbulence  with 
varying  r.m.s.  values  of  intensity  is,  of  course,  still  valid. 
The  intensity  values  of  the  various  Gaussian  patches  are 
characterized  hy  a  probability  density  distribution  q(o) 
with  properties 


cu 

J  q(o)dd  = 


CO 

J  o2q(o)do  =  er^ 


;  (17) 


The  variable 
and  16;  the 


corresponds  to  the  a  used  in  Equations  15 

A 

given  by  Equation  IT  refers  to  the  r.m.s. 


value  of  all  the  patches  linked  together.  With  q  chosen, 
we  use  Equations  15  and  16,  with  the  subscript  x  dropped, 
and  arrive  at  the  composite  results  given  by 


- p(x) 


(18) 


00 

t?(x)  =  PH0  j  q{a)< 
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>'-5^  °  ''' 

C  ^  ?  -f  .*-  r#s  -  - 
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As  is  proper*  it  may  "be  shown  that  the  -composite  density 
curve  mast  satisfy 


cy 

*  J  P'(Jc).dac 

—  CO 
CO 

=  jT  x~p(x)sbc 


P*3r 


Dual  hut  equivalent  definitions  for  the  composite  r.mus. 
value  of  x  are  noted  to  he  given  hy  Equations  17  and  :20„ 


To  illustrate  this  procedure*  we  use  the  well-known 
case  where  q  ^.s  assumed  to  he 


q(a) 


2  1 


v  a 


Substitution  of  this  equation  into  Equation  19  gives  the 
familiar  result 


m  L  i 


U  =  PNQe 


By  summation  timewise ; 

Reference  1  Indicated  that  equivalent  composite  results 
could  also  he  obtained  hy  assuming  that  the  cr^  in  Equation  15 

is  expressed  hy  a  timewise  variation  instead  of  hy  a  distribu¬ 
tion  q  .  In  this  case  the  composite  results  are  given  by  the 
equations 


1  1  Pi 

pW=^'^FI  J  » 

*1 


m  Wm 

fjw  ||3F4$ 

I  pi 
3  %P 
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N(x)  = 


t2“tl 


*&ere  o  is  sene  chosen  function  of  tins, 
As  a  single  example,  suppose 

^=5 


Substitution  into  Equation  22  yields  the  result 


x^tr 


S(x)  ~  nnQ  t^rt^  (e 


2a  ^ 
) 


Tn s  constant-  a  is  established  from  the  fact  that  the. 
value  is  ^luen  by 


.m.s. 


a  '  =  ^rq-p-  / 

x  v2  ±  - 


r  _2 


2  a  J2 

o„  =  ? — —  log  x- 

*  t. 


Additional  Moment  Properties.-  Further  information  on 
moment  properties,  specifickLly  the  relationships  between  the 
various  order  moments  of  the  area  under  the  p  ,  q  ,  and  N 
curves  about  the  origin,  can  be  derived  by  suitable  manipula¬ 
tions  of  Equations  18  and  19 .  Ee cause  of  the  symmetry  of  the 
p  and  if  curves,  odd-order  moments  are  taken  using  the  right- 
hand  plane  onlyj  even-order  moments  consider  both  right  and 
left  halves.  Thus,  from  Equations  .18  and  19,  the  following 
equalities  may  be  derived 


f  sv 
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}  x  pfx)  dx 


far  1 


it  <r  % 

J  1 


f  X*~~;  dx  =  fr  f  xS f)/x/ Jx  -  {Ttt  f  <rJ 


(<r)  dr 
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Analytical  Forms  for  N  and  p  As  derived  from  the 
foregoing  "considerations,,  a  number  of  specific  analytical  ex¬ 
pressions  for  the  "universal 11  curves  for  N(x)  and  p(x)  are 
presented  in  this  section.  Three  forms  are  considered. 

General  form: 

A  large  number  of  specific  analytical  expressions  were 
derived  by  means  of  Equation  10  or  through  use  of  Equations  18 
and  19.  Some  of  the  more  interesting  results  obtained  are 
shown  in  Table  I  as  the  first  nine  entries.  Figure  1  gives 
plots  of  the  associated  exceedance  curves.  For  some  of  the 
eases*  q(o)  is  also  shown  on  the  figures.  The  first  three 
cases  are  noted  to  be  the  same  cases  as  those  given  in 
Reference  4* except  that  here  results  are  given  in  generalized 
form. 


A  comparison  of  the  various  cases  is  given  in  Figure  2; 
the  specific  qurve  chosen  from  each  case  was  arbitrary.  This 
comparison  indicates  that  all  the  curves  may  be  considered 
to  represent  various  deviations  about  a  basic  reference  curve 

defined  by  e  °x  .  If  a  given  universal  curve  falls  above 
this  reference  curve  at  low  x  values*  then  it  must  fall 
below  at  high  x  value  Sj  if  below  at  low  x  *  then  it  must 
be  above  at  high  x  .  In  general*  the  position  of  the  "tail"* 
or  portion  of  the  curve  at  high  values  of  x  ,  is  seen  to  be 
very  sensitive  to  the  initial  position  or  trend.  This  tendency 
is  simply  reflecting  the  property  given  by  Equation  14.  The 
fact  that  the  curves  appear  to  have  large  deviations  from  one 
another  at  large  x  is*  of  course*  due  to  the  use  of  semi¬ 
log  plots j  the  effect*  however*  is  of  much  importance  in  load 
exceedance  and  design  considerations. 

Composite  forms: 

In  the  development  of  the  universal  exceedance  curves 
of  Figure  1  and  in  taking  into  consideration  the  form  suggested 
by  flight  data*  it  was  observed  that  a  desirable  curve  would  be 

X 

-  737 

one  which  fell  essentially  along  the  e  -curve  for  a  wide 

range  of  x  and  which  would  then  depart  in  upward  fashion 
from  this  trend  at  large  x  .  A  universal  curve  of  special 
interest  and  having  these  characteristics  was  therefore 
developed  through  use  of  a  composite  technique.  The  curve  is 
shown  in  Figure  3  and  ,is  derived  from  Equations  13  and  14  as 
follows. 

Assume  £he  upper  and  lower  portions  of  the  curve  are 
given  respectively  by  simple  exponential  and  power  laws  ac¬ 
cording  to  the  expressions 


-ax 

e 


x  <  xi 


b 

y  a  *  x  > 


16 


and  assume  that  the  two  portions  Join  with  equal  height  and 
slope  at  the  joining  point  .  These  two  joining  conditions 
allow  for  the  solution  of  a  xand  b  arid  yield  the  result 


*  x  <  (?.3a) 


\exj 


r-  > 


(2Jo) 


We  now  invoke  relation  14  to  solve  for  x^  5  and  find  that 


(23c) 


Condition  13  is  automatically  satisfied  by  the  exponential 

function  choice.  Equations  23  thus  form  the  universal  curve 

shown  in  Figure  3.  The  equations  and  associated  curves 

applying  for  various  values  of  a  axe  given  below  and  in 

Figure  4.  where  f  and  f.  are  used  in  shorthand  way  to 

a  D  w 

denote  the  two  portions  of  : 

x 

a  =  6, 


a  =  7} 


-1.000818  $ 
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.994274a  n 7  /2.57a  n 
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A  study  of  this  set  of  universal  curves  in  the  light  of  flight 
data  indicates  that  only  the  larger  values  of  a  are  of  prac¬ 
tical  concern.  For  these  larger  values  it  is  noted  that  x, 
is  ^iven  in  good  approximation  by  ~ 


occr„ 


and  thus  Equations  23a  arid  23b  read  simply 


(24a) 


N  _  faoA 


,  ~  >  “  (24b) 

X 


Strictly  speaking,  to  preserve  the  moment  relation  given  by 
Equation  14,  the  coefficient  of  unity  in  the  exponent  of 
Equation  24a  should  be  minutely  larger  than  one  (as  given  by 
230)5  the  difference  between  plotted  results  when  using  unity 
or  the  precise  value  cannot  be  detected,  however,  and  so  the 
use  Ox'  unity  for  practical  purposes  is  justified.  The  combi¬ 
nation  of  Equations  24  thus  forms  a  remarkably  simple  and 
easy  to  apply  universal  curve  for  load  exceedances.  Its  ease 
of  application ' will  be  seen  in  some  of  the  subsequent  design 
applications.  Equations  24  also  include  another  case  that  has 
been  of  special  interest;  that  is,  when  cc  =  «>  the  result  is 
simply 
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This  is  noted  to  be  Case  a  of  Table  1. 

Cases  k  and  i  represent  two  other  interesting  ccrpo- 
site  forms  that  were  derived  through  ass  of  Sqaaticns  13  2nd 
14  and  which  show  characteristics  similar  to  Case  3  Just 
discussed;  plots  are  given  in  Figures  3  and  6.  Case  k  is 
noted  to  be  a  refined  version  *>f  the  fos 


'  ?2e 


which  has  been  referred  to  often  in  previous  studies.  The 
newer  developments  of  the  present  report  have  been  used  to 
express  the  equation  in  a  form  that  is  compatible  with  the 
other  cases. 

Superposed  forms: 

After  performing  the  work  covered  in  the  preceding 
sections,  which  had.  the  underlying  aim  of  developing  reasonably 
simply  analytical  expressions  for  the  exceedance  curves,  and 
which  at  the  same  time  yielded  desired  shape  characteristics, 
the  following  fact  was  realized.  Any  number  of  different  ex¬ 
ceedance  curves  can  be  added  together  by  coefficients  which 
sum  to  unity  to  yield  yet  another  exceedance  curve.  Speci¬ 
fically,  in  terms  of  two  different  chosen  exceedance  curves 
f-y  and  ■>  ■fcke  following  applies  in  general 


Case  m  of  Table  I,  and  associated .Figure  7>  illustrate  this 
technique.  This  case  was  formed  from  Cases  a  and  d  and  is 
of  special  interest.  It  Is  quite  simple  analytically,  and, 
in  the  light  of  available  flight  exceedance  data.,  shows  good 
shape  characteristics . 

With  the  variety  of  cases  presented  in  Table  I,  the 
question  naturally  arises  as  to  shat  forms  are  roost  appro¬ 
priate  for  practical  use.  In  partial  answer  to  this  question, 
we  must  admit  that  at  the  present  time  not  enough  experimental 
data  exist  in  true  generalized  exceedance  form  to  allow  a  se¬ 
lection  to  be  made  (or  to  associate  form  with  mission  type). 
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&.  study  sabseqaazrtl;;  in  this  report  using  a.  liEited  asoam; 

of  flight  data  indicates  tthich  of  the  forzss  appear  to  he  fa¬ 
vored.  On  the  oasis  of  the  inforaaticoi  available  at  the  pre¬ 
sent  tise.  Cases  J  and  e  ^re  preferred  choiCne. 


\'Tr* 
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SECTION  4 

BSJ5IC  PARAMETERS  FGR  DESIGN 


In  general  there  are  five  basic  parameters  which  are 
involved  in  the  consideration  of  gust  loads.  Two  are  the 
basic  structural  response  parameters  A  and  2?Q  *  and  three 

are  the  basic  environmental  parameters  brought  out  by  Equa¬ 
tion  11 2  namely  P  ,  °x  ~  Aoe  5  a  *  ^4ese  parameters, 

with  the  exception  of  a  ,  are  discussed  In  this  section. 
Because  of  their  close  tie-in,  A  and  a  are  discussed 

together.  A  sixth  parameter,  the  scale  of  turbulence  L  , 
may  also  be  mentioned.  The  scale  L  and  the  parameters  A 
and  dc  are  so  interrelated,  however,  that  only  two  of  the 

three  may  be  considered  independent.  Thus  comments  on  the 
scale  value  will  also  be  made  when  appropriate  in  the  dis¬ 
cussion  of  A  and  .  A  further  discussion  of  scale,  par¬ 
ticularly  the  effect  that  filtering  of  time  history  gust  data 
has  on  deduced  scale  value,  is  given  in  the  appendix. 

All  of  the  subsequent  discussion  is  based  on  the  fact 
-  which  is  demonstrated  hy  Equation  11  that  each  one  of  the 
parameters  NQ  P  ,  and  ax  may  be  considered  to  be  a  com¬ 
posite  value  of  $11  the  flight  experience.  These  composite 
values  may  be  established  conveniently  by  assuming  that  the 
flight  of  an  aircraft  is  expressed  in  terms  of  various  flight 
segments  as  depicted  in  the  following  sketch: 


Segment  1 
dci 


Segment  2 
9c  p 


Segment  3 


The  segments  are  arbitrary  and  may  represent,  for  example, 
climb,  cruise,  and  descent  -  or  any  other  convenient  and  ap¬ 
propriate  breakdown.  The  model,  somewhat  analogous  to  the 
discrete-gust  patch  model  concept  discussed  in  Reference  3* 
corresponds  to  a  condensed  version  of  a  mission  approach.  In 
the  development  of  this  report,  the  determination  of  exceedance 
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curves  for  eacL  segment  is  obviated*  yet  generality  and  a  feel 
for  the  sensitivity  of  mission  results  to  the  individual  seg¬ 
ments  can  still  be  retained*  as  will  be  seen. 

General  relationships  which  apply  to  this  model  and  which 
are  to  be  used  throughout  the  developments  to  folio*  are 


35* 


0.2  = 


piTl  4  p2^2  +  P3T3 

' 

(^} 

i  (Pl~lacJ  , 4  ?2T2ae^  + 

P3T3a03) 

(^) 

+  + 

E3V*3> 

(27) 

where,  f pr  cogyehieiice , ih.  writing,  consideration  has  been 
.ii^ted:'t^:thr%e  ~s'eg^nts, r~  -  - 

„  Additional  basic  relationships  involving  the  structural 
p^rametdr r  A  ' '  "  '  are  *  '  •  •  : 


.  i 


(28b) 


The  Frequency  Parameter  The  total  number  of 

upward~iero  crossings  for  the  model  would  be 


n0  P2T2N02  +  P3T3N03 


From  this  equation  the  composite  follows  simply  as 

N0  -  ff  ■  M  SlTl\  +  Wo2  +  P3T3H03)  W 

The  basic  spectral  formula  for  evaluating  the  individual 
N0*s  is  given  in  Reference  1. 
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The  Composite  Prooortica  of  Tims  Parameter 


past  efforts  save  attempted  to  establish  the  average  amount 
of  time  that  -is  spent  in  turbulence,  and  particularly  the 
variation  with  altitude .  So  suits  vary  considerably,  Referen¬ 
ces  §,  6,  and  7-  Because  different  analytical  techniques 
were  used,  there  is  scene  question  oh  the  reliability  of  the 
results.  In  an  attempt  to  ascertain  appropriate  values  of 

P  -  wsa  mans  TOri  fnie  sms  —  -SI  <->+• o 


experience  involving  crop  dusting,  charter,  and  corporate 
flying,  presented,  a  picture  oh  proportion  of  time  which  is  con 
sidered  truly  amazing.  He  verified  that  geographic  location 
and  season,  were,  of  course,  dominant  factors  influencing  gust 
encounter.  For  an  overall  coverage,  however,  Mr.  Murphy 
suggested  a  breakdown  involving  daytime  (8  a.m.  -  8  p.m.)  and 
nighttime  (8  p.m.  -  8  a.m.),  operation.  The  results  he  gave, 
purely  from  mental  reflection,  are  shown  in  Figure  8.  Since 
these  restilts  are  not  dependent  on  a  specific  analytical  re¬ 
duction  procedure,  they  are  considered  the  most  reliable  yet 
shown.  '  •  ' 


As  an  illustration  of  how  the  values  of  P  shown  in 
Figure  8  may  he  used  In  mission  consideration,  consider  the 
daytime  values  to  be  represented  by  the  solid  straight  line 
segments  which  are  defined  as  follows; 


T?  =  ;8  e 


-  .1623  - 


1000 


-  .0886  1QQ0 


0  <  h  <  15,000 


15,000  <  h  <  30,000 


30,000  <  h 


These  equations  may  be  combined  with  any  assumed  mission  in 
accordance  ’with  Equation  25  to  yield  the  effective  value  of  P 
Because  of  the  analytical  form  assumed  for  the  .P's  ,  exact 
integration  can  be  used  in  place  of  a  summation  convention. 
F^r  example,  the  results  found  from  Equations  30  and  the  as¬ 
sumed  mission 


\ 

a 

T  1 

are  the  following 


0  <  h  <15*000,* 
p  “’8  5a  + 


•8’  1  "  Tra.  O  +  ***)] 


15,000  <  h  <  30,000, 


l1  -  (  1  +  15,000a)] 


30,000  <  1), 


p  -  -8  &  *  & .  ( I  -  1  +  !&«»*)  + 


•1  -  Jb  (1  +  wo} 


and  where 


is  the  rate  of 


where  a  =  *  b  =  and  where  r  is  the  rate  of 

climb  and  descent*  here  considered  equal.  Figure  9  shows  the 
results  obtained  from  these  equations.  These  results  repre¬ 
sent  expected  average  values  of  P  for  the  assumed  mission 
for  daytime  operation.  Because  of  the  similarity  of  day  and 
nighttime  curves  shown  in  Figure  8*  estimates  of  p  for  opera¬ 
tions  which  represent  a  combination  of  day  and  night-time 
flights  may  he  made  quite  easily  from  Figure  9  through  use  of 
a  simple  reduction  factor.  Thus*  for  a  combination  of  half 
daytime  and  half  nighttime  flights*  P  may  be  estimated 
simply  by  multiplying  the  'ralue  shown  in  Figure  9  by  .6  . 

For  all  nighttime  operation  a  factor  of  .3  appears  reasonably 
good. 

The  opinion  is  given  here  that  values  of  P  for  opera¬ 
tions  which  use  turbulence-avoidance  procedures,  as  practiced 
generally  by  the  commercial  airlines*  are  probably  not  too 
much  less  than  the  values  given  in  the  preceding  paragraph. 
Usually  the  intent  of  avoidance  procedures  is  to  avoid  mainly 
the  severe  turbulence*  and  this  type  contributes  only  a  small 
amount  to  P  .  The  principal  effect  of  avoiding  the  severe 
turbulence  probably  shows  up  more  as  a  change  in  the  value  of 
the  shape  parameter  a  ,  as  examples  later  in  the  report  will 
show. 

Results  for  P  were  given  in  Reference  1  for  a  different 
assumed  basic  variation  with  altitude*  but  no  equations  were 
given.  The  variation  considered  was  given  by  the  exponential 
form 
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p  =  pco 


Consideration  of  the  basic  elimb-eruise-deseent  type  mission 
led  to  the  following  equation: 


(*  -  >)  *  0-  -  S)  •'“* 


where  Ir  is  the  mission  cruise  altitude,  ibis  equation 

gives  the  results  shown  in  Figure  4  of  Reference  1,  Vol.  II 
(the  subscript  m  is  dropped  on  both  P  and  h  in  the 
figure).  The  results  given  are  for  the  assumed  values  of 
PQ  =  .2  and  a  «  .08/1000  . 


The  Composite  Severity  Parameter 


and  the  Related 


Parameter  A  .-  The  composite  severity  parameter  cr  if 

—  A 

given  hy  Equation  27.  From  the  basic  relations  28,  this 
severity  value  becomes 


a 2  = 
x 


w(wf°0i +  +  W3ac  ZJ -  a2°2 


In  application,  the  A  !s  and  cr  's  of  this  equation  cannot 

n  cn 

he  treated  independently  at  the  present  time  because  they  are 
closely  linked  together  by  the  turbulence  scale  value  L  . 

The  values  of  A  are  strongly  dependent  on  the  scale  value 
used  in  the  chosen  gust  spectral  equation  arid,,  in  the  con¬ 
sideration  of  flight  data,  it  has  not  been  possible  thus  far, 
because  of  instrumentation  limitations,  to  separate  the  ef¬ 
fects  of  a  and  L  *  Thus,  neither  a  or  L  ,  or  the 

variation  of  these  quantities  with  altitude,  are  really  known 
at  this  time.  %e  must  await  better  flight  test  data  to  as¬ 
certain  their  values.  In  the  meantime,  the  approach  given  in 
Reference  1,  making  use  of  the  truncated  value  of  gust  inten¬ 
sity,  is  suggested. 

The  truncated  value  of  intensity,  as  associated  with  the 
area  tinder  the  spectrum  beyond  a  frequency  ^  }  is  given  by 

(Rote,  ct  and  ii&  here  replace  the  and  as  used  in 


Reference  1,  and  g.„  Is  the  gust-  intensity  of  one  of  the 

ft 

patoises  that  make  up  the  composite  value  .)  In  general, 

is  given  with  good  accuracy  by  trie  flight  data,  and  so 

the  hast  that  can  be  said  about  cl.  and  L  is  that  they 

** 

are  seme  cccsbinsticn  that  satisfies  this  relation .  Without 
sore  additional  piece  of  Information,  no  unique  values  can 
be  established,  (Soes  bounds  exist,  of  course,  on  o._  %  i.e., 

R 

it  cannot  be  smaller  than  o*.  nor  Is  It  larger  than  Hie  mean 

square  value  of  the  rasr  time  history  data  of  gust  velocity.) 
Lacking  specific  information  on  a.  ,  we  can  treat  q*.  as  the 

basic  gust  intensity  variable.  Concepts  of  a  distribution 
curve  for  ,  similar  to  that  used  in  Equations  18  and  lg, 

and  a  variation  of  e-  with  altitude,  similar  to  that  given 

xc 

for  o  in  Figure  3  of  Beferenee  1,  Vol.  H,  are  equally  valid; 
Os.  represents  the  eomncsite  value  of  oi  ,  analogous  to  the 

v  ~  r  -  -  —  i> 

C  . 

use  of  c  as  the  composite  value  of  0  .  We  may  write, 

therefore. 


.335<jc 
(LJV)^/3  * 

v 


(33) 


Chi  the  basis  of  the  limited  information  that  exists,  and  also 
as  partially  guided  by  the  results  that  were  obtained  in  the 
reexamination  of  some  airline  data  as  given  later  in  this 
report,  we  assume  the  following:  that  0.  =  2,45  for  an 

zc 

associated  =  o003  •  These  assumptions  lead  to  the  g& 

values  shown  in  Figure  10.  This  figure  is  tentatively  sug¬ 
gested  for  design  purposes. 


An 

lowing.  If 
we  find  that 


equivalent  procedure  for  treating  oQ  is  the  fol- 


is  eliminated  between  Equations  28b  and  33, 
>1/3 


0X  =  A 


.885 


Vt 


(34) 


Reference  1  indicates  that  is  fairly  insensitive  to  L  . 

Thus,  as  an  alternative  to  the  use  of  Figure  10,  design  couli 
make  use  of  Equation  34,  with  o ^  =2.45  fox  =  .003  as 

the  tentative  choice  for  the  truncated  gust  intensity  vslue. 
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In  this  form  the  necessity  of  considering  L  is  virtually 
eliminated . 


The  composite  value  of  A  is  of  interest  and  may  be 
derived  directly  .from  Equations  2 6  and  32;  the  result-  is 


.2  _  ^L.  =  ^iAi  *  Vf  +  ^3 

~  2  7i  +  7o  +  7o 


(35a) 


where  = 
'n 


*A‘  2 


If  a  =  cr  =*  a  ,  as  is  sug- 
C1  c2  c3 


gested  in. Reference  1,  and  in  the  preceding  paragraphs,  the 
solution  for  A  is  simply 


,s_*A 


,2  +  a2  + 

j.  ^  PT '  "2  PT 


(35b) 


Thus,  Equations  35  provide  for  the  fairly  easy  evaluation  of 
the  effective  or  composite  value  of  A  for  use  In  load  ex¬ 
ceedance  considerations.  An  illustration  of  the  application 
of- Equation  35b  to  a  fairly  complex  type  mission  is  shown  in 
Figure  11;  the  numbers  given  near  the  bottom  of  the  figure 
show  which  of  the  segments  contribute  most  to  the  composite 
value  of  A  , 

With  respect  to  the  concept  of  using  an  increased  a, 

w 

value  when  increased  severe  turbulence  encounter  is  contem¬ 
plated,  it  should  be  mentioned  that  the  results  given  in 
Figure  5  of  Reference  1,  Vol.  II,  are  slightly  in  error;  the 
correct  values  should  be  the  square  root  of  the  values  shown. 
The  equation  appropriate  to  the  figure  is  derived  as  follows. 
Consider  the  sequential  encounter  of  turbulence  of  intensity 
for  a  time  T^  and  of  a2  for  a  time  Tg  ;  the  total  time 

of  turbulence  encounter  is  thus  T  =  +  T2  .  The  composite 

mean  square  value  is  given  by 

4  - 

“fC/Wiat  +  /w2at) 

0  T-T2 


-  rji  (TjC-j.  +  T2C2) 

1 
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The  correct  values  of  ac/cr^  as  obtained  from  this  equation 
are  given  in, Figure  12  (note  and  Pp  of  Figure  5  in 

Reference  1,  Vol.  II,  are  replaced  here  respectively  by  T^ 
and  Tp) . 

Errors. and  Sensitivity  to  A  and  The  sensitivity 

of  N  “to  possible  errors  in  A  or  Nn  ,  or  to  errors  in  cr 

v  v 

and  P  ,  is  of  interest.  One  way  to  show  this  sensitivity 
is  to  assume  that  A  or  NQ  may  be  in  error  by  a  given  per¬ 
centage  and  to  determine  how  much  N  has  changed  as  a  result. 
Figure  13  shows  typical  results;  the  figure  shows  the  percentage 
change  in  N  due  to  an  assumed  10  percent  error  in  either  A 
or  oc  ,  or  Nq  or  P  .  In  Figure  14,  error  effects  are  brought 

out  in  another  way;  number  of  exceedances  as  might  be  obtained 
in  -a  specific  case  are  shown  for  base  values  of  A  ,  ,  oc  , 

and  P  ,  and  for  values  of  these  parameters  which  are  assumed 
to  be  in  error  t  10  percent.  The  main  effect  noted  from  these 
figures  is  that  small  errors  in  A  or  a  can  lead  to  large 

errors  in  N  .  A  10  percent  error  in  either  or  P  gives 

a  10  percent  error  in  N  ,  but  a  10  percent  error  in  A  or 
a  can  lead  to  errors  in  N  of  several  hundred  percent. 

These  observations  emphasize  the  fact  that  A  and  oc  are  the 

critical  parameters  and  need  to  be  established  as  accurately 
as  possible,  but  that  in  comparison,  approximate  estimates 
suffice  for  and  P  .  It  is  noted  that  errors  in  A  or 

a  can  be  likened  to  the  use  of  a  different  a  curve.  Thus 
c 

the  choice  of  a  and  the  accuracy  by  which  A  can  be  estab¬ 
lished  are  closely  interrelated. 


SECTION  5 

EFFECT  OF  MISSION  OH  SHAPE  OF  EXCEEDANCE  CURVE 


The  four  possible  design  procedures  that  were  advanced  in 
reference  1  are  to  he  discussed  further  in  a  subsequent  section 
on  DESIGN,  The  subject  of  mission  design  merits  some  additional 
clarifying  thoughts  however,  and  this  consideration  is  given  in 
this  section.  The  mission  design  procedure  that  is  being 
followed  here,  as  in  reference  1,  represents  a  departure  from 
that  explored  in  the  past.  Here  we  do  not  necessarily  build  up 
the  mission  curve  by  adding  together  the  curves  from  all  the 
various  assumed  flight  segments.  We  simply  choose  a  single 
universal  curve  to  represent  the  mission.  For  example,  in 
Figure  4,  the  curve  for  a  =  8  may  be  stipulated  as  the  curve 
applying  to  Mission  A  ,  while  the  curve  a  =  6  may  be  chosen  to 
represent  Mission  B  .  Thus,  in  the  mission  approach  under  con¬ 
sideration  here  we  simply  specify  one  of  the  uni;  al  curves 
(which  actually  may  be  established  by  experience  ._i  another 
aircraft)..  In  addition,,  we  need  values  of  c>c  and  P  .  These 

may  be  established  in  two  ways:  (l)  by  stipulation*,  or  (2)  they 
may  be  calculated  through  mission  segment  considerations  by 
Equations  25  and  26  (note  a  and  P  may  be  established  on  a 

V 

segmented  basis  even  though  the  N  curves  are  not  handled  in 
this  way).  With  the  universal  curve,  and  P  chosen,  the 

design  curve  Is  completely  specified  (except  for  the  value  of 
turbulent  scale  L  ,  which  is  handled  as  discussed  in  connection 
with  Equations  33  and  34. 

As  a  matter  of  Interest,  the  effect  of  mission  on  the  shape 
of  the  mission  universal  curve  ia  demonstrated  by  the  following 
treatment o  Consider  a  mission  composed  basically  of  three 
different  segments,,  an  ascent,  cruise,  and  descent  portion.  We 
may  construct  a  universal  curve  applying  to-  the  whole  mission  in 
a  manner  analogous  to  the  segmented  mission  approach,  but  here, 
vje  assume  that  each  segment  is  governed  by  the  same  basic  uni¬ 
versal  curve,,  say  one  of  the  curves  in  Figures  4  or  7*  We  write 
for  convenience  that  N(  ) /PNq  =  f(  );  the  composite  curve  for’ 

load  exceedance  then  follows  as 


n  -  plTlN01f(o~)  +  +  p3 wG“) 


(37) 


X2' 


The  analysis  and  equations  that  apply  for  determining  the  com¬ 
posite  values  P  ,  a  ,  a,.  ,  and  N.^  are  given  in  the 

v  A  \ 

preceding  section.  Through  these  equations  Equation  37  may  be 
written  in  the  universal  form 
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(38) 


where  ^  =  (^V^^OiA)} 


Altitude  Effect.  -  To  see  how  the  shape  of  the  universal 
curve  might  be  influenced  by  mission.  Equation  38  was  applied 
to  the  assumed  high  altitude  mission  shown  in  the  following 
sketch  s 


Note  that  the  gust  velocity  is  assumed  to  be  invariant  with 
altitude.  Results  obtained  using  the  a  =  8  curve  of  Figure  4 
as  the  basic  universal  curve  are  shown  in  Figure  15.  The  a  =  8 
universal'  durve  is  shown  also  and  is  labeled  zero  altitude,  since 
this  curve  is  precisely  the  result  that  would  be  obtained  for  an 
assumed  "zero  altitude™  mission.  A  comparison  of  results  and 
referral  to  Figure  4  indicates  that  the  main  effect  of  altitude 
is  to  yield  a  universal  curve  having  a  slightly  lower  value  of  a 

Turbulenc e  Intensity  Ef f ec t .  -  The  effect  on  shape  of 
Increased “operation in " severe  type  turbulence  can  be  demonstrated 
in  a  similar  way.  Consider,  for  example,  an  operation  having 
turbulent  encounter  as  shown  in  the  following  sketch 


c. 


If  this  model  is  analyzed  in  a  manner  analogous  to  that  used 
for  Equation  38,  the  following  result  is  obtained 


r/*jv\ .  !a  y 

[A*®,,  V  Ti  V 


'*  V 
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(39) 


where 


is  given  by  Equation  36,  and  where 


fc  fl 
°1  g2 


Figure  16(a)  shows  the  universal  curves  that  are  obtained 
in  this  case  for  various  values  of  Tg/T^  and  0^/0 ^  , 

arbitrarily  choosing  the  Case  (a)  universal  curve  of  Figure  1 
as  the  basic  curve.  With  reference  to  Figure  4,  we  see  that  the 
tendency  to  shift  to  curves  having  lower  values  of  a  is  much 
more  pronounced  in  this  instance  than  for  the  altitude  case. 
These  tendencies  and  trends  suggest  that  a  reasonable  scheme  for 
selecting  the  universal  curve  to  use  in  design  application  could 
be  of  the  nature  indicated  by  the  following  listing: 

Case  i  Case  m 

1.  Normal  low  altitude  mission  a  =  8  .005 

2.  Normal  high  altitude  mission  a  =  7.5  ,01 

3.  Extended  flight  through  severe 

turbulence  a  =  6  .04 

The  numbers  given  here  are  not  necessarily  recommended  choices; 
they  are  intended  mainly  to  indicate  the  general  trends  or 
behavior.  It  should  be  noted  that  the  discussion  here  Is 
centered  mainly  on  the  shape  of  the  curve;  each  mission  may  have 


different  composite  values  of 


and  they,  of  course 


individually  have  a  marked  influence  on  the  load  exceedance 
history,  as  much  or  more  so  than  the  shape  factor.  In  the 
examples  given,  these  values  were  suppressed  in  the  role  of 
normalizing  factors  so  as  to  emphasize  the  shape  effect. 

Figure  16(b)  is  an  extension  of  the  results  of  Figure  16(a) 
to  show  that  as  the  assumed  operation  approaches  complete  severe 
turbulence  encounter,  T.,  0  ,  the  universal  curve  returns  to 

the  same  form  as  that  for  no  severe  turbulence  encounter, 

Tg  0  .  The  load  history  for  the  high  ratios  of  T2/T^  will, 

of  course,  be  much  more  severe  than  for  the  low  ratios  of  Tr./T1 


because  <?c  is  greater,  even  though  the  shapes  of  the  universal 

curves  are  essentially  the  same  as  for  the  case  =  0  * 

pie  most  marked  change  in  shape  is  noted  to  occur  srtien  only  a 
small-  amount  of  severe  turbulence  is  involved . 
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SECTION  6 

HBEXflMIHATIQN  OF  OPERATIONAL  DATA 


The  applicability  of  the  universal  curves  of  Figure  4 
was  tested  by  applying  the  curves  to  the  flight  operational 
data,  given  in  Reference  4.  Part  of  this  reexamination  con¬ 
sisted  of  reevaluating  the  values  of  A  and  that  applied 

to  the  airplane s,  because  only  erude  estimates  of  these  para¬ 
meters  were  used  in  the  reference  report 5  the  rigid  body  re¬ 
sults  given  in  Section  5  of  Reference  1,  Vol.  II,  were  used 
for  this  purpose.  Results  of  the  test  are  shown  in  Figure  17; 
associated  values  of  the  various  parameters  that  were  evalua¬ 
ted  and  deduced,  such  as  P,  0,0,  are  listed  in  Table  II. 

The  fitting  of  the  universal  curves  given  by  Equation 
24  and  pictured  in  Figure  4  involved  two  kinds  of  fit:  (1) 
at  least  two  of  the  leftmost  data  points  fell  on  the  straight 
line  section  of  the  universal  curve.  Operations  1,  2,  3>  4, 

5,  6;  (2)  only  one  point  fell  on  the  straight  portion  so 

that  the  Joining  point  of  the  two  portions  was  between  the 
first  and  second  point.  Operations  7  and  8. 

The  results  for  A  and  ac  in  Table  II  (a)  apply  to 

an  assumed  scale  L  =  1000  ft.  In  Table  11(b)  the  results 
for  other  choice  of  L  are  given  as  well.  The  wide  varia¬ 
tion  shown  in  the  deduced  value  of  a  ,  depending  on  the 

value  of  L  used,  demonstrates  the  role  and  importance  of 
knowing  the  scale  value.  As  another  interesting  point  to  be 
seen,  the  table  shows  also  the  influence  of  flight  route  on 
the  amount  of  time  spent  in  turbulence;  specifically,  the 
largest  values  of  P  are  seen  to  be  associated  with  geogra¬ 
phic  locations  where  more  turbulence  is  generally  to  be 
expected. 

Some  comments  on  the  success  of  the  reexamination 
follow.  The  matching  of  the  flight  data  with  the  universal 
curve  chosen  is  perhaps  the  best  ever  seen.  There  is  no 
arbitrariness  about  fitting  the  curves.  Further,  it  appears 
that  the  curves  provide  a  consistent  and  reliable  means  for 
deducing  the  values  of  P  and  a  that  apply.  It  is  noted 
also  that  in  the  cas^  of  Equations  24,  the  properties  of 
exponential  and  powers  laws  are  such  as  to  make  the  fitting 
of  the  analytical  expressions  to  the  data  a  rather  easy  task. 
With  the  encouraging  success  indicated  here,  it  is  considered 
desirable  to  examine  in  future  work  the  applicability  of  the 
various  generalized  exceedance  curves  to  other  flight  data. 

Many  aids  may  be  devised  to  simplify  the  task  of  fitting 
a  chosen  universal  curve  to  flight  data.  One  is  to  plot  ex¬ 
ceedance  values  against  x^  on  log-log  paper  to  Judge  whether 
the  data  tend  to  behave  as  a  powei  law  at  high  x  ;  the  use  of 
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the  fourth  power  is  arbitrary,  but  the  .choice  helps  to  reduce 
the  steepness  of  the  curves  and  aids  in  fairing.  Another  ai- 
is  the  following:  convert  the  exceedance  values  to  a  N/80 

basis,  divide  the  variable  x  by  several  assumed  values  of 

ax  ,  multiply  the  W/N0  values  by  e_X//<7x  and  plot  the  results 

against  x/c?  .  The  value  of  a  which  causes  the  data 
points  to  lie  along  a  hor  i  z ont alxl ine  at  low  x  is  the  de¬ 
duced  flight  value;  the  value  where  this  horizontal  line  in¬ 
tersects  the  vertical  axis  is  the  P  valued  To  find  the  a 
value,  multiply  the  chosen  universal  curves  also  by 

g-x/O-x  and  plot;  then  lay  the  flight  results  over  these 
modified  analytical  curves  so  that  the  chosen  flight  hori-  ■ 
zontai  line  overlays  the  horizontal  line  for  the  analytical 
curves  (which  should  be  unity) .  The  value  of  a  is  then 
found  by  interpolation  by  noting  where  the  flight  values  at 
high  x/o_  fall  relative  to  the  various  analytical  curves. 

X  % 

This  reexamination,  using  the  curves  of  Figure  4,  was 
made  before  the  curves  of  Figure  7  (and  some  of  the  other 
cases)  had  been  established.  It  is  expected,  however,  that 
equal  success  and  similar  P  and  a  values  would  have  been 
found  if  the  curves  of  Figure  7  had  been  used. 

A  summary  of  the  generalized  exceedance  curves  as  de¬ 
rived  from  the  flight  data  is  given  in  Figure  18  as  a  matter 
of  interest.  We  conclude  this  section  with  an  observation 
that  is  in  the  nature  of  a  warning .  We  should  not  draw  any 
hard  conclusions  from  this  reanalysis,  or  from  the  summary 
given  in  Figure  18  for  the  following  reasons.  First,  only 
a  limited  amount  of  data  has  been  considered  here.  Secondly, 
the  exceedance  data  used  is,  strictly  speaking,  not  of  the 
proper  kind;  the  data  represents  a  certain  type  of  peak 
count  analysis,  but  whether  the  specific  count  analysis  used 
represents  the  results  that  would  be  obtained  by  a  level 
crossing  analysis  -  on  which  the  exceedance  curves  are  based  - 
is  not  known.  Therefore,  we  cannot  judge  at  this  time  which 
family  (or  families)  of  generalized  curves  represents  flight 
exceedance  the  best.  We  should  also  keep  in  mind  the  following 
thought  in  any  future  analysis:  any  exceedance  analysis  should 
also  include  data  at  low  x  values  because  knowledge  of  the 
behavior  of  the  exceedance  data  at;  low  and  intermediate  x 
can  be  very  valuable  in  predicting  how  the  exceedance  curves 
should  behave  at  very  large  x  ,  where  data  is  usually  non¬ 
existent.  Stated  in  another  way,  the  behavior  of  the  ex¬ 
ceedance  curve  at  large  x  is  very  strongly  dependent  on 
whether  the  curve  at  low  values  of  x  is  convex,  straight, 
or  concave,  when  viewed  from  above  (on  the  basis  of  a  semi-log 
plot) . 


SECTION  7 

rssiGs  modssm&s 


In  the  light  of  the  newer  developxoents  that  have  been 
found  in  this  phase  of  study,  we  examine,  in-  this  section, 
the  four  design  procedures  that  were  advanced  in  Reference 
Essentially,  the  four  procedures  still  appear  sound. 

Some  of  the  numbers  are  altered  slightly  herein;  other  num¬ 
bers  still  cannot  be  fixed  with  certainty  until  better  and 
more  appropriately  analyzed  flight  data  become  available. 

'She  discussion  is  presented  for  the  most  part  using  a 
general  notation  for  uhe  universal  exceedance  curve.  In 
places.  Equation  24  is  used  to  give  specific  illustration. 

In  the  lifetime  T  of  the  airplane,  the  number  n  of 
exceedances  of  level  x  is  (see  Equation  11) 

n  =  TN(x)  =  PTH0f  (4 
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where 


Most  of  the  discussion  to  follow  is  given  on  the  basis  of 
limit  load  considerations;  a  treatment  of  ultimate  load  Is 
also  given,  however,  tc  show  that  some  insight  on  the  near¬ 
ness  to  or  possibility  of  ultimate  load  encounter  can  be 
gained,  even  though  possible  complications,  such  as  nonlinear 
structural  and  aerodynamic  behavior,  are  not  considered.' 

At  limit  load,  x  =  Ax^  ,  we  wish  the  expected  number 

of  exceedances,  as  given  by  Equation  40,  to  be  equal  to  or 
less  than  a  specified  number  n^  ;  thus 


In  terms  of  Equation  24b,  this  relation  reads  specifically 

^(sr)  <  "jb 


35 


(42) 


Design  Based  on 


*0 


vs.  x/A  For  safe  design.  Equa¬ 


tions  40  and  41  indicate  the  following  functional  relation 
between  NQ  and  .x/A  : 


*0< 


_ n  t/x\ 

*  “) 


With. -E quat ion  40,  this  would  read  specifically 


Bo£°l 


<£ 


(*3) 


As  mentioned  in  Reference 


•values  to  n 


T 


1,  c-,  can  be  established  by  assigning 

,  and-1  a  ,  or  it  may  be  found  in  em¬ 


pirical  manner  by  applying  the  equation  in  retrospect  manner 
to  past  gust-critical  airplanes  that  have  proved  themselves 
airworthy  by  many  years  of  luccessful  operation.  Let  us  exa¬ 
mine  what  the  second  of  these  two  ways  indicates.  In  Figure 
1 9,  we  plot  the  results  of  the  computational  studies  that  are 
presented  in  References  8,  9,  and  10.  We  see  the  exoected 
^shot  gun"  pattern  of  points  and  note  that  they  establish  a 
reasonably  well-defined  left-hand  border,  as  was  anticipated 
and  hoped.  Equation  43  was  used  to  establish  the  three  border 
curves  shown;  the  curves,  representing  a  =  7>  8,  and  9,  were 
all  arbitrarily  made  to  pass  through  the  point  x/A  =  52.5  3 
Nq  =  1  .  Specifically,  the  chosen  borders  are  given  by  the 

equation 


N, 


-  te)° 


(44) 


It  is  noted  that,  because  of  the  general  steepness  of  the  curves, 
the  border  itself  is  not  too  critically  dependent  on  the  value 
of  a  chosen.  Note  also  that  quite  similar  borders  won.1!  be 
established  if  other  specific  forms  of  the  exceedance  curve 
had  been  used. 


The  results  shown  in  the  figure  are  for  computations 
based  on  a  scale  L  =  1000  ft.  If,  however, the  results  for 

x  / 1000 

other  L  values  are  plotted  using  (  -j- — )  as  the 

abscissa,  as  based  on  the  indications  of  Equation  34,  then 
similar  patterns  with  the  same  left-hand  boundary  are  found. 

It  would  appear  desirable,  therefore,  to  use  the  abscissa 
which' includes  the  L  factor,  since  the  necessity  for  specifying 
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L  thereby  becomes  obviated. 

An  important  point  to  note  In  this  procedure  is  the 
following*  In  the  consideration  of  various  flight  conditions* 
it  is  possible  that  one  or  more  flight  conditions  may  give 
combinations  of  SfQ  and  x/A  for  a  specific  structural  de¬ 
sign  point  that  falls  to  the  left  or  unsafe  side  of  the  de¬ 
sign  border*  while  all  other  flight  conditions  yie3  d  values 
which  fall  to  the  right  or  safe  side.  The  composite  value 
combination  of  and  x/A  *  as  given  by  Equations  29  and 

35b  may*  however*  still  fall  on  the  safe  side.  A  procedure 
as  follows  is  therefore  suggested.  Evaluate  and  x/A 

for  various  flight  conditions*  and*  if  the  points  all  fall 
on  the  safe  side*  the  design  may  be  fudged  safe;  there  is  no 
need  to  establish  composite  or  mission  values.  If  some  of 
the  points  fall  on  the  unsafe  side*  then  determine  the  compo¬ 
site  values  of  and  x/A  ;  in  all  likelihood*  these  will 

fall  on  the  safe  side*  thus  indicating  a  safe  design.  If  the 
composite  values  fall  on  the  left  side*  the  design  should*  of 
course*  be  judged  unsafe. 

Figure  20  presents  some  of  the  airplane  computational 
results  to  illustrate  these  situations.  For  the  RB-57*  stress 
point  9  appears  safe  for  all  four  flight  conditions  assumed. 
For  stress  point  13  however,  a  combination  for  and  x/A 

is  obtained  which  falls  just  to  the  left  of  the  desi^i  border. 
The  flight  condition  for  this  case  is  5$  fue3  at  50,000  ft. 
altitude.  In  this  situation  the  design  may  be  judged  safe 
without  establishing  the  composite  values  of  and  x/A  * 

since  the  amount  of  time  spent  in  this  apparent  critical 
flight  condition  is  surely  quite  small.  A  similar  situation 
is  found  to  occur  for  stress  point  1  for  the  B-58*  whereas 
for  stress  point  3a  results  fall  on  the  safe-  side  for  both 
flight  conditions  investigated.  The  results  shown  for  the 
KE-135  ar®  seen  to  be  safe*  and  serve  also  to  show  the  spread 
in  points  due  to  various  flight  conditions. 

In  general*  the  concept  of  using  past  aircraft  to  es¬ 
tablish  a  design  border  for  this  vs.  x/A  approach  ap¬ 

pears  successful.  The  use  of  composite  values  of  NQ  and 

x/A  when  necessary  makes  the  approach  quite  attractive  and 
versatile.  The  nature  of  the  results  found  suggest*  in  fact* 
a  simplification  worthy  of  consideration.  As  noted*  the  de¬ 
sign  borderline  is  characterized  by  a  marked  steepness*  so 
much  so  that  it  appears  a  vertical  line  design  border  might 
serve  just  as  well  for  practical  purposes.  From  the  results 
shown  in  Figure  19*  we  might  choose*  for  example*  a  vertical 
line  at  an  abscissa  value  of  55.  Design  could  then  be  couched 
in  the  following  concise  phrasing;  design  for  gusts  shall  be 
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sues  as  zo  satisfy  -r  I 

course,  the  tj  c  approach  a c 


>  55  .  This  approach  is,  of 


coarse,  the  tio„  autrrc-acb  advanced  in  Reference  3l  it  is 

analogous  to  tbe  discrete-gust  approach  and  baa  simplicity 
as  one  of  its  merits .  These  results  also  bring  out  again  a 
point  that  was  discussed  In  a  previous  section,  namely,  that 
A  is  predcginant  in  design,  and  that  relatively,  has  a 

secondary  role. 

Design  Based  on  Comparison.-  The  comparison  chart  that 
follows  iron  Equation  40  and  which  takes  the  place  of  Figure 
6  in  Reference  1,  Vol.  H,  is  given  in  Figure  21.  The  curves 
shown  are  based  on  Equation  42.  This  method  is  still  con¬ 
sidered  good.  Hot  only  is  it  good  for  comparing  one  aircraft 
design  to  another  design,  dug  it  should  he  especially  good 
for  comparing  the  adequacy  of  the  design  at  one  point  of  an 
airplane  with  the  design  at  another  point  of  the"  same  airplane . 
This  latter  point  deserves  much  attention.  Thus,  some  other 
procedure  may  he  used  to  evaluate  or  fix  the  design  at  a 
given  ei'tical  reference  point  on  the  wing,  for  example,  and 
then  use  can  he  made  of  the  comparison  technique  to  checx. 
other  points  along  the  wing  or  elsewhere .  When  applied  in 
this  way  to  the  same  design,  there  is  no  concern  over  what  the 


values  of  P  ,  T  ,  and 


(or  B  )  should  he.  Again,  com¬ 


posite  values  as  given  hy  Equations  25,  2 6,  29,  and  35b,  may 
he  used  if  desired  or  needed. 

Design  Based  on  Stipulated  <jc  ,  P  ,  and  T  In  this 
procedure,  we  specify  the  values  of  cj  ,  P  ,  T  ,  and  L  to 

v» 

he  used  for  design  (and  a  particular,  universal  exceedance 
curve) j  we  also  fix  the  number  n^  of  limit  load  encounters 

acceptable.  Values  of  A  and  NQ  are  computed  for  various 

flight  conditions  and,  for  simplicity,  may  first-  he  checked 
individually.  If  some  are  found  to  he  on  the  unsafe  side  for 
certain  flight  conditions,  then,  as  in  the  previous  two  ap¬ 
proaches,  composite  values  of  A  and  NQ  (and  P  )  may  he 

used  to  test  acceptability  from  a  mission  point  of  view.  De¬ 
sign  can  he  examined  completely  and  directly  on  the  chosen 
universal  exceedance  curve,  as  is  indicated  in  Figure  22. 

For  limit  load.  Equation  4l  indicates 
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The  left-hand  side  of  this  relation  fixes  the  level  of  the 
horizontal  line  a  .  For  a  design  to  he  considered  satis¬ 
factory,  level  a  must  he  equal  to  or  above  point  A  . 
Ultimate  load  encounter  may  also  he  treated,  at  least  in  a 
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roi  gh  sense,  by  introducing  the  number  m  of  airplanes  in  the 
fleet,  We  do  not.  wish  to  encounter  an  ultimate  load  during 
the  lifetime  of  the  fleet  and  thus,  on  the  basis  that  Equation 
40  also  applies  in  approximate  fashion  to  ultimate  loads,  we 
have 

nu  -  1  >  -«V  (§»«). 


In  this  case,  the  left-hand  side  fixes  the  level  cf  the  hori¬ 
zontal  line  b  .  The  example  numbers  shown  in  Figure  22  serve 
to  illustrate  the  procedure. 

The  “master  design"  chart  form  suggested  in  Reference  1 
may  also  be  used  in  an  alternative  way  to  assess  the  design. 
Figure  23  illustrates  the  form,  using  Satiation  24b  with 
a  =  8  to  establish  the  curves.  This  particular  form  was 
suggested  because  of  the  analogy  to  and  possible  use  for 
fatigue  considerations,  as  brought  out  in  Reference  1. 

The  following  comments  are  intended  to  be  of  a  general 
nature,  and  indicate  one  means  for  establishing  the  number 
of  limit  load  encounters  that  might  be  appropriate  for  design 
use.  The  process  is  admittedly  open  to  some  question  .since 
ultimate  load  aspects  are  involved,  and  since  the  numbers 
found  are  dependent  on  the  particular  exceedance  curve  chosen., 
but  at  least  order  of  magnitude  values  should  be  indicated. 

If  Equation  24b  and  the  abscissa  relations  shown  on  Figure 
22  are  used,  the  following  simple  relation  for  the  ratio  of 
the  number  of  limit  load  encounters  to  the  number  of  ultimate 
load  encounters  is  found 


n 


u 


This  relation  yields  the  results 
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Statistically,  these  mashers  indicate  the  number  of  limit  load 
encounters  that  may  be  expected  for  the  fleet  for  one  expected 
ultimate  load  encounter ;  for  example,  512  are  indicated  for 
the  case  of  r  =■=  1  and  a  =  9  *  -Prom  these  numbers  and  a 
consideration  of  airplane  fleet  size,  we  can  establish  the 
number  of  limit  load  encounters  per  airplane  that  is  reasonable 
to  use  for  design.  Thus,  for  the  example  ratio  of  512,  the 
following  relation  between  fleet  size  m  and  limit  load  en¬ 
counters  n^  per  airplane  must  be  preserved  so  as  not  to  ex¬ 
ceed  (statistically)  one  ultimate  load  encounter  of  the  fleet 
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These  numbers  singly  reflect  the  fact  that  if  more  airplanes 
are  involved,  then  the  chances  of  encountering  ultimate  load 
increase .  Other  exceedance  curves  would  indicate  different 
numbers,  but  we  must  await  better  experimental  evidence  to 
judge  which  curve  is  the  most  realistic  to  use. 

Design  Based  on  Mission  Approach  (Utilizing  exceedance 
curves  for  each  aeg^nt).-  The  mission  approach,  in  which  ex¬ 
ceedance  curves  fS:aeit£sS)lished  for  each  assumed  mission  seg¬ 
ment  and  then  aumj&ed,  is  not  developed  further  herein  as  a 
design  tool.  A  suppl amenta!  report  to  this  report  is  inten¬ 
ded,  however,  fee  cover  additional  developments  and  newer 
ideas  on  this  approach .  In  reading  the  subsequent  discussion, 
the  following  thought  should  be  kept  in  mind.  Most  of  the  ob¬ 
servations  are  made  hare  in  a  negative  vsi n;  positive  aspects 
will  be  dealt  with  in  the  supplemental  report.  The  lack  of 
support  .Indicated  bore  is  not  because  the  mission  approach  is 
felt  to  be  improper.  The  principal  investigator  of  the  present 
report,  in  fact,  first  advanced  the  use  of  generalized  ex- 


ceea&nce  curve  in  a  mission  approach  in  Reference  11,  and 
restated  the  concept  in  Reference  3-  The  approach  is  still 
considered  good,  hut  unfortunately  not  enough  information 
exists  as  yet  to  place  it  on  an  acceptable  firm  basis. 

At  the  present  time  the  segmented  mission  approach  is 
recommended  not  as  a  direct  design  tool,  but  rather  as  a 
means  for  assessing  how  the  character  of  the  expected  load 
history  is  changed  if  the  aircraft  is  utilised  in  different 
ways.  ’There  are  several  reasons  for  these  observations. 

The  main  reason  is  that  too  many  uncertainties  still  exist 
with  respect  to  the  input  quantities.  In  the  treatment  of 
physical"  problems  in  which  -uncertainty  exists,  there  is  some¬ 
times  a  tendency  to  add  yet  another  ingredient  to  the  problem 
in  the  belief  that, the  situation  is  made  clearer.  In  reality 
though,  the  problem  has  been  made  all  the  more  uncertain. 

At  present,  the  segmented  approach,  while  potentially  a  very 
powerful  tool,  has  many  of  these  aspects.  In  particular, 
when  evaluated  on  the  basis  of  the  number  of  statistical 
variables  involved,  the  approach  has  four  independent  variables 
these  may  be  stated  in  various  ways,  but  a  common  way  expresses 
them  as  2  proportion  of  time  values,  and  2  gust  severity 
values.  While  attempts  have  been  made  to  establish  these 
values,  they  are  still  very  uncertain. 

The  approaches  given  herein  can  accomplish  everything 
that  can  be  done  with  the  segmented  approach,  but  in  a  much 
abbreviated  way,  and  with  fewer  variables.  Corresponding  to 
the  four  variables  mentioned,  there  are  only  three,  or  less, 
involved  in  the  approaches  of  this  report  -  specifically  P  , 
ac  ,  and  a  (for  the  vs.  x/A  approach,  there  are  none). 

The  uncertainty  of  a  fifth  variable,  which  is  common  to  all 
procedures,  the  scale  L  ,  has  also  been  removed  to  a  large 
extent.  In  a  sense  the  uncertainties  that  exist  with  all  the 
specific  generalized  exceedance  curves  of  the  mission  approach 
are  here  lumped  together  into  a  single  curve  described  by  the 
parameter  a  .  The  underlying  thought  is  that  instead  of 
summing  together,  say,  10  different  curves,  each  of  which  has 
a  degree  of  uncertainty,  why  not  choose  a  single  curve  repre¬ 
senting  the  sum,  which,  even  though  it  too  may  be  somewhat 
uncertain,  is  at  least  guided  or  constrained  more  directly 
by  past  experience  results. 

An  illustration  of  the  pitfall  that  may  be  encountered 
if  design  is  made  on  a  mission  basis  is  afforded  by  back 
reference  to  Figure  lfe.  These  results  were  obtained  from 
three  different  airplanes  of  the  same  type  that  were  utilized 
in  three  different  ways,  but  at  essentially  the  same  mission 
altitude.  The  question  suggested  by  these  results  1st  how 
mi{  it  a  mission  have  been  conceived  for  design  that  could 
have  covered  the  load  histories  actually  experienced?  A 
mission  covering  the  lowest  curve  might  have  been  conceived 
only  to  find  that  the  results  given  by  the  top  curve  were 
obtained  in  utilization  of  the  aircraft.  The  answer  is  that 
the  mission  conceived  must  cover  all  possible  future  uses  of 
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SECTION  8 

CONSIDERATION  OP  STRUCTURAL-  INTERACTION  EFFECTS 


The  exceedance  considerations  discussed  so  far  in  this 
report  apply  to  any  response  quantity  which  can  be  treated 
separately.  It  is  known,  however,  that  there  are  often  points 
in  aircraft  structures  where  the  strength  is  governed  by 
structural  interaction  boundaries;  that  is,  the  strength  de¬ 
pends  on  the  combination  of  two  stresses,  such  as  an  axial 
stress  and  a  shear  stress."  The  problem  of  determining  the 
number  of  times  a  given  interaction  boundary  is  exceeded  is 
therefore  present. 

The  means  for  determining  the  e,xpected  number  of  ex¬ 
ceedances  of  a  given  structural  interaction  boundary  has  been 
developed  in  a  separate  study.  The  procedure  Is  general  but 
quite  simple,  and  exact  solutions  have  been  obtained  for  cer¬ 
tain  cases.  A  report  covering  the  study  is  under  preparation 
and  thus  no  further  discussion  will  be  given  here. 
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SECTION  9 

PROBABILITY  OF  EXCEEDANCE  FOR  GIVEN  FLIGHTS 


In  the  development  of  this  report  there  is  no  explicit 
consideration  given  to  establishing  probability  of  failure 
or  the  probability  levels  of  exceeding  certain  loads,  such  as 
limit  load.  The  reason  is  that  there  is  no  need  to  do  this, 
There  may  be  interest  and  concern,  however,  in  the  question 
of  what  is  the  probability  of  exceeding  a  certain  load  3.evel 
in  a  given  number  of  hours  of  routine  flight,  or  for  flights 
such  as  extended  operation  through  known  severe  turbulence 
conditions.  An  analysis  of  this  probability  problem  is  given 
in  this  section.  -The  results  are  basic  and  apply  in  general 
to  the  probability  of  failure  or  the  probability  of  survival 
for  structures  exposed  „o  random  loading.  Three  solutions 
are  presented. 

A  fundamental  ingredient  to  all  solutions  is  the  average 
repeat  time  for  a  given  load  level.  From  the  generalised  ex¬ 
ceedance  Equation  11,  this  value  is 

m  _  1  _  1 

x  “  N  “  PNQf  (45) 

jf  x  \ 

where  f  is  understood  to  represent  the  function  ff  — a  }  . 

\cx  / 

We  ask  then,  what  is  the  probability  that  the  load  level  x 
will  be  exceeded  in  time  T  . 

The  first  solution  is  approximate.  We  assume  that  T 
is  very  small  compared  to  Tx  .  Then  the  probability  of  ex¬ 
ceedance  is  simply 

.  -  P(x,T)  =  S- 

-V  xr 


=  TPN0f  (46) 

The  solution  obviously  breaks  down  when  T  approaches  or 
exceeds  Tx  .  (Note  in  this  preser.tat.ion  the  distinction 
between  the  functional,  notation  P(x,T)  for  the  probability 
and  the  pro-portion  of  time  P  should  be  kept  in  mind;  P  wil 
be  used  throughout  she  remainder  of  this  section  to  denote  the 
proportion  of  time  P  .) 

The.  second  solution  follows  the  line  of  reasoning  used 
in  a  private  communication  from  Prof.  B.  Etkin.  Let  T  be 
divided  into  a  sequence  of  small  intervals  At  such  that 


With  Equation  45  this  yields  the  desired  result  for  P 


-  TPNQf(x} 


P(x,T)  =  1  -  e 


This  equation  is  noted  to  he  associated  with  the  Poisson 
distribution. 

It  is  questionable  whether  Equation  47  properly  repre 
sents  the  actual  characteristics  of  a  random  function.  For 
one  thing  the  equation  is  based  on  the  average  repeat  time, 
but  in  reality  the  repeat  time  is  variable  and  must  be 
characterized  by  a  distribution  function.  Secondly,  the 
equation  does  not  yield  the  correct  results  when  applied' 
to  a  deterministic  function  such  as  a  sine  wave,  as  is  the 
belief  that  ic  should.  A  development  which  overcomes  both 
of  these  objections  is  the  following.  Consider  the  random 
function  shown  in  the  following  sketch 


nSa 


We  are  interested  in  the  probability  of  whether  the  level  x 
will  be  experienced  in  a  given  time  interval  T  .  By  way  of 
illustration  we  first  analyze  the  situation  on  the  basis  that 
t  is  characterized  by  a  discrete  value  distribution;  after¬ 
wards  a  continuous  distribution  function  for  t  iss  considered. 

A  *  ''**■  ■ 


Suppose  that  the  individual 


t  values  are  read  over  a 

A 


long  time  interval  and  that  when  these  values  are  grouped  In 
interval  brackets,  that  the  following  specific  distribution 
diagram  is  found 


tl  t2  t3  t4  t5 


To  find  the  probability  of  encountering  x  in  time  T  we 
can  make  use  of  a  probability  wheel  of  the  following  type 


The  factor  a  ±u  inserted  simply  to  convert  the  time  values 
to  angle  notation.,  such  that 


i 


Assume,  for  illustration  purposes,  that 

t3  <  T  <  t^ 


Then,  through  the  consideration  of  many  repeated  spins  of  the 
dial,  the  probability  that  the  shaded  area  will  not  cover  an 
x  encounter  in  a  single  try  is  • 


9fx  T)  -  g(t4  -  T)  +  (t5  -  T) 

t  +  J5t  +  3t0  4  2t  i.  4  t_ 

1  2  ~  3  4  5 


(48) 


The  probability  of  encountering  x  is  thus 

2(t4  -  T)  +  (t  -  T) 
P(x,T)  -  1  Q(x,T)  1  “  ti  +  2tg  +  3t,  4  2t4  4 
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The  average  repeat  time  is  found  from  the  first  moment  of 
the  xi  diagram,  or  specifically 


t±  +  2tg  -r  3t3  +  2t^  -5-  t5 
Tx  =  1  +  2  +  3  +'  S'+T 


(^9) 


if 


*  iAT 


A  generalization  of  this  elementary  model  may  he  made 
as  shorn  hy  the  following  density  plot  and  equations 


P 


where 


n  =  ^  +  ng  +  n^  +  . . , 


Si  - 1 

i 


(50) 


In  analogy  to  Equations  48  and  49>  and  for  T  <  ,  we  may 

write  the  probability  of  not  encountering  x  as 


Q(x,T) 


ni  ^i“T^  +  ni+l 


T)  +  n-*+P  (^+2~^)  + 


i+2  '  i+2 


n^t.  +  n0tg  +  n^E^T* 
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and  the  average  repeat  time  T ,  as 


fp 

~x 


+ 


05 


'fi 

3-1 


> 


t 


(5‘ 


In  the  limiting  case  as  the  interval  AT  =  becomes 

small,  the  p  distribution  may  be  expressed  in  continuous  form 
as  shown  in  the  following  sketch 


which  by  Equations  50  and  52  has  the  properties 


The  generalized  Q  equation  becomes 

oo 

/  ft  -  T)p(x,  t)at 

Q(x,T)  =  - - - -  =  7jr;  f  (t  -  T)p(x,t)dt 

J'  tp(x,t)dt  ^ 

0 

We  have  finally,  therefore,  that  the  probability  of  encountering 
x  in  a  time  T  is  given  by  the  equation 


P(x,T)  -  1  -  Q(x,T) 


=  1 


J  (t  -  T)p(x,t)dt 
T 


(54) 


This  relation  is  general  and  should  apply  to  all  stationary 
functions  whether  random  or  not.  The  precise  form  of  p  de¬ 
pends  on  the  nature  of  the  random  variable ;  for  a  sine  wave 
p  becomes  a  Dirac  function  at  the  period  of  the  wave , 

The  nature  of  p  for  gust  encounter  loads  is  not 
known  as  yet,  but  estimates  of  P  can  nevertheless  be  made 
by  assuming  various  forms  of  p  ,  Two  choices  have  been 
studied j  they  are 

Case  a,  _  nt 

=  Txr(n)  (f~)  '  e  * 


Case  b, 

p(x,t)  =  ate 


The  second  case  is  presented  in  form  only  because  the  coeffi¬ 
cients  a  and  b  are  defined  by  rather  lengthy  expressions 
involving  exponentials  and  erf(  ).  Figure  24  shows  the 
graphical  nature  of  these  distribution  functions.  The  results 
found  from  Equation  54  for  these  two  cases  are  shown  in  Figure 
25,  along  with  the  results  indicated  by  Equations  46  and  4?, 
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To  bring  out  better  the  behavior  of  the  curves  near  P  =  1  * 
probability  of  not  exceeding  are  given  in  Figure  26.  The 
results  obtained  for  extreme  values  of  either  n  or  k  are 
noted  to  be  of  special  and  remarkable  interest.  Case  a,  yields 
Equation  46  for  n  =  «  ,  and  Equation  47  for  n  =  1  ,  whereas 
Case  b*  gives  Equation  46  for  k  =  <»  and  for  t ,  =  -<»  yields 
the  n  -  2  result  for  Case  a.  x 

.  The  results  are  general.  In  application  to  the  gust 
loads*  it  is  noted  that  the  results  are  presented  in  terms 
of  the  average  repeat  time  as  given  by  Equation  45. 

Although  Information  is  not  available  at  this  time  to 
judge  the  validity*  the  results  given  by  Equation  54  look 
very  satisfying.  The  associated  curves  given  in  Figure  25 
appear  as  generalizations  to  Equation  47.  It  Is  noted  also 
that  the  distribution  functions  p  approach  a  Dirac  function 
as  the  parameters  n  and  k  become  large*  and  hence  apply 
to  a  sinusoid  as  a  limiting  case;  the  corresponding  proba¬ 
bility  curves  pass  in  the  limit  to  the  approximation  given 
by  Equation  46*  which  Is  correct  for  a  sine  wave ,  Although 
there  is  no  supporting  evidence,  it  is  felt  that  the  proba¬ 
bility  curves  in  the  range  of  n  =  10  -  20  In  Figure  24(a) 
and  k  =  2  -  2.5  in  Figure  24(b)  apply  to  the  gust  load  ex¬ 
ceedance  case . 

The  results  given  in  Figure  25  can  be  used  in  several 
ways  to  determine  the  probability  of  exceeding  certain  loads. 
One  of  course*  is  on  the  basis  of  routine  flying.  Perhaps  a 
more  significant  use*  however*  is  with  respect  to  establishing 
probabilities  for  a  certain  flight*  or  succession  of  flights, 
into  known  severe  turbulence  conditions.  Results  can  be  quite 
startling.  Two  examples  illustrate  the  point.  Consider  that 
a  design  was  made  using  the  following  values: 

Ax. 

*■* j ”  16 

P  =  .08 

T  =  15*000  hrs,  =  54  x  106  sec. 

N0  =  1.05 

«L  =  5 

and  that  the  a  ~  .005  curve  of  Figure  7  was  the  universal 
exceedance  curve  used  for  design.  The  k  =  2  curve  of 
Figure  25(b)  is  chosen  to  determine  probabilities. 

Example  1: 

We  wish  to  establish  the  probability  of  encountering 


a  load  equal  to  half  the  limit  load  In  10  hours  of  routine 
flight 5  for  this  instance  we  have 


t(jf)  ~  f  (8)  ~  4  x  10" 


§-  «  PTN0f^-|^)  =  .  08  r  10  x  3600  x  1,05  x  4  x  10"^ 


=  1.21 


From  Figure  25(h)  we  find  P  j=  .955.  In  addition  we  wish 
now  to  determine  the  time  to  reach  this  same  probability 
level  if  flight  Is  made  completely  In  turbulence  that  Is 
2,5  times  as  severe  as  the  average.  In  the  latter  Instance 

P  =  1,  ox  is  2.5  times  as  large  or  =  3.2,  f (3.2) 

1  -2  T  x 

-  4.2  x  10  ,  and  sp-  must  be  the  same  to  yield  the  same 

A 

probability,  hence 


PTN0f  ( )  =  1.05T  x  4.2  x  10"2  *  1.21 


or 


T  =  27.5  sec. 


In  summary,  for  this  example.  It  takes  10  hours  of  normal 
routine  flying  to  yield  a  probability  of  ,955  that  a  load 
level  one-half  the  limit  load  value  will  be  reached,  but  only 
27.5  sec.  of  flying  In  turbulence  which  Is  2.5  times  as  se¬ 
vere  as  the  average  to  achieve  the  same  probability. 

Example  2: 

We  wish  to  determine  how  much  time  can  be  flown  con¬ 
tinuously  in  turbulence  that  Is  2.5  as  severe  as  the  average 
to  obtain  the  same  probability  of  encountering  limit  load 
as  is  obtained  in  flying  routinely  for  a  time  equal  to  the 
design  life  of  the  airplane.  For  the  lifetime  routine  flight 
we  have 


~~  =  PTNrm6)  :=  .08  x  54  x  106  x  1.05  x  1..  x  Hf6 
x 

~  b 


The  associated  probability  is  practically  1,  For  the  severe 
turbulence  flight  to  have  the  same  probability  we  have  *  with 

T  =  1,  — ~  =  tti?  -  6.4,f(6.4)  =  1,66?  x  10~3,  the  following 
x 

relation 


5  =  PTH0f 


~  1.05T  x  1,667  x  10~3 


and  so 

T  -  47.7  min. 


These  examples  serve  to  show  the  value  of  a  probability 
of  exceedance  analysis,  and  indicate  that  probability  consider a 
tlons  of  the  type  Illustrated  can  be  quite  useful  in  Judging 
whether  a  specific  flight  in  known  severe  conditions  should 
be  made  or  not. 


SECTION  Xv 

CONCLUSIONS  AND  RECOMMENDATIONS 


The  conclusion?  and  r  commendations  that  are  indicated  by 
this  study  dealing  with  the  development  of  power  spectral 
methods  for  designing  aircraft  for  gusts  are  listed  in  this 
section. 

First,  as  a  general  observation,  the  use  of  spectral  tech¬ 
niques  for  gust  design  appears  soimd’and  fe asjf^le,  and  it  is 
recommended  that  su"h  a  procedure  be  Incorporated  in  design 
specifications.  This  recommendation  does  not  mean  that  the 
discrete-gust  design  technique  should  be  abandoned.  The  atti¬ 
tude  taken  toward  the  discrete-gust  method  is  that  because  of 
its  basic  simplicity  it  still  may  be  used  to  rough  out  the 
design  in  the  early  stages;  the  power  spectral  technique  then 
serves  fchs  purpose  of  checking  on  the  design  or  uncovering 
unusual  response  aspects  which  cannot  be  brought  out  rationally 
by  the  discrete-gust  technique.  We  might  summarise  the  use  of 
the  power  spectral  approach  as  follows;  (l)  If  sufficient 
structural  and  aerodynamic  information  is  available  early  in  the 
design  stage,  then  design  by  the  power  spectral  approach,  (2) 

If  the  phasing  of  the  design  is  such  as  to  preclude  making  use 
of  the  spectral  approach  initially,  then  rough  out  the  design 
by  the  discrete-gust  procedure  and  subsequently  make  a  detailed 
check  by  the  power  spectral  approach,  (3)  E vert  if  not  used 
initially  for  a  given  design,  the  power  spectral  method  repre¬ 
sents  a  powerful  tool  1  j  determine  how  the  expected  load  history 
may  be  influenced  or  changed  if  the  airplane  is  used  in  different 
ways  or  in  a  manner  different  from  that  originally  contemplated. 

With  respect  -co  future  study  effort,  the  following  observa¬ 
tions  are  made;  7 

1)  The  values  of  a  ani  L  are  still  uncertain  and 
should  be  established  more  rationally. 

2)  The  questionable  reliability  of  the  derived  gust 
history  records  at  low  frequencies  is  one  of  the 
principal  reasons  for  the  state  indicated  inlfcem 
u)i  thus,  a  much  better  appreciation  of  record 
content  must  be  established. 

3)  The  operational  statistics  of  various  aircraft 
should  be  examined  to  evaluate  which  form  (or 
forms)  of  the  many  generalized  exceedance  curves 
presented  appears  to  represent  load  experience 
the  best. 

4)  The  study  mentioned  in  item  3  should  also  Include 
ar  examination  of  whether  the  shape  of  the 
exceedance  curves  can  be  correlated  with  mission 

category. 
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5)  Load  level  crossing  studies  should  he-  made  of 
gust  response  records-  results  should  he  cor¬ 
related  with  previous. j.,;-  used  peak  count  proce¬ 
dures,  Any  level  crossing  stuuy  should  also 
include  the  range  of  low  x  valuer-,  since  in¬ 
formation  in  this  range  may  provide  a  valuable 
clue  as  to  how  the  tall  of  the  exceedance  curves 
should  behave. 

6)  Additional  values  of  and  x/A  should  be  es¬ 
tablished  for  other  gust  critical,  aircraft  as  a 
means  for  providing  further  substantiation  to  the 
results  and  design  boundary  given  in  Figure  19* 

7}  Gust  response  records  should  be  evaluated  to  es¬ 
tablish  experimental  distribution  functions  for 
repeat  time  as  used  in  the  probability  study 
given  in  Section  9» 


With  respect  to  the  spectral  design,  procedures;, 
following  recommendations  are  made* 

1)  The  design  based  on  the  plot  of  Nn  vs. 


the 


1000 


1/3 


’o 

and  represented  by  Figure  19*  is  the  recommended 
first  choice  as  a  design  procedure.  Composite 
mission  values  of  A  and  N0  may  be  used  if 

necessary.  This  procedure  is  just  about  as  simple 
as  can  be  conceived,  yet  which  retains  all  the 
basic  rational  elements  of  the  problem.  Serious 
consideration  should  be  given  to  adopting  the 
further  simplification  of  this  procedure  -  that 

of  specifying  the  lower  limit  to 

only,  which  on  Figure  19  Is  a  vertical  line. 


2)  Design  based  on  the  universal  curve  and  specified 
(or  composite)  values  of  oQ  ,  P  ,  and  T  (Figures 

22  or  2g)-  is  recommended  as  a  second  choice.  As 
mentioned,  uncertainties  still  exist  as  to  what 
values  of  on  ,  P  ,  and  L  should  be  used.  The 

v 

procedure,  however,  provides  for  additional  rele¬ 
vant  design  information,  such  as  load  Information 
for  fatigue  considerations. 


3)  The  comparison  technique  is  not  included  in  an 
order  of  choice  basis.  This  procedure  is  consi¬ 
dered  a  valuable  adjunct  to  any  procedure  that  Is 
being  used.  It  is  valuable  in  evaluating  the  de¬ 
sign  at  various  locations  on  a  given  aircraft, 
and  in  many  instances  can  itself  serve  to  be  the 
complete  design  tool. 
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4)  The  segmented  mission  approach  involving  the  use 
of  an  exceedance  curve  for  each  mission  segment 
is  not  ia -eluded  as  a  design  procedure  at  this 
time.  The  main  reason  is  that  too  many  uncer¬ 
tainties  still  exist  with  respect  to  the  input 
Quantities.  The  use  of  this  procedure  as  a  means 
for  assessing  how  the  expected  load  history  is 
affected  hy  different  aircraft  utilization  is, 
however,  recommended.  Further  developments  on 
this  approach  are  to  do  covered  in  a  supplemental 
report , 


APPENDIX  A 


EFFECTS  OF  FILTERING  ON  L 


Evidence  Is  that  the  time  histories  of  gust  velocities 
as  derived  from  airplane  flight  data  contain  long  wavelength 
contaminations  of  uncertain  magnitude  and  frequency  range. 
Knowledge  of  this  low  frequency  “noise”  is  important  because 
its  presence  has  a  major  influence  on  the  values  of  a  and 

v 

L  that  are  deduced.  As  an  example,  in  Reference  8,  it  is 
shown  that  scale  values  of  over  50,000  ft.  are  found  in  some 
cases  from  the  raw  time  history  records  of  derived  gust  velo¬ 
cities;  such  values  appear  to  be  unreasonably  out  of  line. 
Because  of  the  dependence  of  a  and  L  on  record  authenti¬ 
city,  the  principal  investigator  of  this  report  has  recommen¬ 
ded  often  that  long  still  air  flight  tests  which  include  de¬ 
liberate  pilot  Induced  maneuvers  be  conducted  with  the  instru¬ 
mented  airplane.  The  records  obtained  from  such  flights 
should  give  a  good  indication  of  what  can  be  believed  and  what 
cannot  in  the  gust  records.  In  the  absence  of  these  studies, 
other  means  for  ;cmbating  the  apparent  low  frequency  noise 
have  been  used.  One  of  these  methods  consists  of  arbitrarily 
filtering  out  tie  low  frequency  components.  The  effect  that 
such  filtering  has  on  the  deduced  values  of  scale  should, 
however,  be  appreciated  more.  An  analysis  is  therefor*-  given 
in  the  following,  which  Indicates  the  general  effect  of  fil¬ 
tering. 

Consider  a  raw  time  history  of  a  random  variable  y(t). 
The  time  history  obtained  by  passing  this  raw  function  through 
a  low  pass  filter  is 


00 

yf(t)  =  fy(x)  h(t-r)  dr 
—  00 


where  h  describes  the  filter  characteristics  (although  not 
necessary,  we  assume  that  the  filter  used  has  symmetrical 
characteristics  so  that  pnase  distortion  is  not  introduced) » 
The  time  history  with  low  frequency  components  removed  would 
thus  appear 


00 

I(t)  =  y(t)  -  J  y(r)  h(t  — t)  dr 
**00 


The  Fourier  integral  transform  of  this  equation  is 


=  Py(®)  -  PyH  H(<d) 


where  H  is  the  filter  frequency  response.  The  spectrum  for 
Y  then  follows  as 


*y(«)  1 3  -  H(<o)j2*y 


<A1) 


To  establish  the  effect  of  filtering  on  the  da  cueed 
value  of  scale  L  ,  we  mahe  use  of  Equation  6  of  Reference  1, 
Vol.  II.  giving  the  truncated  r.m.s.  value  .  namely. 


„-885<I 


V 


(&%) 


T/3 


(A2) 


and  postulate  the  equivalence  shown  in  the  following  sketches. 


The  spectrum  on  the  left  applies  to  the  raw  aata,  and  on  the 
assumption  that  these  data  represent  real  gust  information,  we 
fit  our  chosen  analytical  spectrum  equation  (say.  Equation  4, 
Vol.  II,  Reference  1)  and  deduce  an  apparent  or  "raw"  scale  L  , 
The  second  sketch  represents  the  filter  transfer  function.  In 
accordance  with  Equation  Al,  curve  a  in  the  third  sketch 
represents  the  spectrum  of  the  filtered  time  history.  The 

p 

quantity  denotes  the  ^ean  square  value  of  the  filtered 

time  history.  We  assume,  however,  that  spectrum  a  is  re¬ 
placed  by  an  "equivalent"  spectrum  b  ,  on  the  basis  that  b 
represents  the  real  world  situation  much  better  than  a  ; 
stated  from  another  viewpoint,  if  we  had  a  means  for  precisely 
extracting  the  contamination  from  the  raw  data,  then  the  fil¬ 
tered  data  would  presumably  lead  to  spectrum  b  ,  not  a  . 
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The  Fourier  integral  transform  of  this  equation  is 


FyH  =  Py(®)  -  Py(^)  H(a>) 


where  H  is  the  filter  frequency  response,  The  plectrum  for 
Y  then  follows  as 

0y(co)  *-  |.l  -  H(o>)J20y  <A1) 


To  establish  the  effect  of  filtering  on  the  da  aieed 
value  of  scale  L  .  we  mzt he  use  of  Equation  6  of  Reference  1, 
Vol.  II.  giving  the  truncated  r.m.s.  value  oy  ,  namely. 


,885crte 

(i^F3 


(A2) 


and  postulate  the  equivalence  shown  in  the  following  sketches. 


Seale  h 


The  spectrum  on  the  left  applies  to  the  raw  a-ata,  and  on  the 
assumption  that  these  data  represent  real  gust  information,  we 
fit  our  chosen  analytical  spectrum  equation  (say.  Equation  4, 
Vol,  II,  Reference  1)  and  deduce  an  apparent  or  "raw”  scale  L  . 
The  second  sketch  represents  the  filter  transfer  function.  In 
accordance  witn  Equation  Al,  curve  a  in  the  third  sketch 
represents  the  spectrum  of  the  filtered  time  history.  The 

p 

quantity  denotes  the  vean  square  value  of  the  filtered 

time  history.  We  assume,  however,  that  spectrum  a  is  re¬ 
placed  by  an  "equivalent"  spectrum  b  ,  on  the  basis  that  b 
represents  the  real  world  situation  much  better  than  a  ; 
stated  from  another  viewpoint,  if  we  had  a  means  for  precisely 
extracting  the  contamination  from  the  raw  data,  then  the  fil¬ 
tered  data  would  presumably  lead  to  spectrum  b  ,  not  a  , 
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The  equivalent  spectrum  is  defined  as  one  which  has  the  shape 
given  by  the  analytical  expression  adopted*  which  has  a  mean 

p 

square  value  cr®  and  which  passes  through  the  same  data 

points  at  the  high  frequency  end*,  and  hence  has  the  same  trun¬ 
cated  value  '  The  scale  value  Lg  associated  with  this 

equivalent  spectrum  must  then  satisfy  an  equation  similar  to 
Equation  (A2 j *  specifically 


(A3) 


Division  of  Equation  A2  by  Equation  A3  leads  then  to  the  f ol 
lowing  re .ii  ion  defining  the  scale  L„  which  effectively 
chacterizer  the  filtered  data  J 


This  simple  equation  can  be  used  to  show  the  effect  of  filtering 
on  the  deduced  scale  value.  Results  are  shown  in  Figure  27  for 
two  different  assumed  ideal  filters.  The  following  equations 
define  the  filters  chosen 

.64 3o>c  -.643a>c|t| 

=  2'~  e 


Case  1 


h 


Case  2 


h  = 


1  sin 
IT  t 


1  -  H  (co)  j 2  =  0 


O)  <  CO, 
CO  >  CO 

V 


n 


and 


/I* 

Note,  ft  =  —  ,  ftc  corresponds  to  o>c 


0) 

OK 


Note  that  the-  scale  L-  associated  with  the  filtered 
data  drops  off  rapidly  as  the  filter  cutoff  frequency  £1  is 

y 

increased,  and  more  significantly  that  at  the  higher  ft.  1 b  the 

V 

values  of  Lf  become  the  same  regardless  of.  the  value  of  the 

raw  scale  L  .  There  is  perhaps  a  significant  implication  in 
these  results.  If  we  consider,  for  example,  that  the  informa¬ 
tion  contained  in  the  raw  data  for  frequencies  below  ft  =  .0005 

(wavelengths  greater  than  X  =  Sr/il  =  12,560)  is  all  contain!- 

y  y 

nation,  then  we  see  that  scale  values  on  the  order  of  1200  ft. 
are  indicated.  We  see .from  these  results  the  importance  of 
establishing  whether  the  information  given  by  the  raw  time 
history  is  real;  or  not.  The  re  stilts  also  suggest,  the  warning: 
do  not  arbitrarily  filter  the  data  unless  definite  guidelines 
are  on  hand  to  guide  the  filtering  and  unless  the  effects  of 
filtering  on  the  end  results  are  fully  understood. 
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Figure  25b.  Probability  of  Exceeding  Load  Level  x  in  Time  T  5  Case 
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Figure  26.  Probability  of  Not  Exceeding  Load  Level  x  in  Time 
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established.  It  is  shown  that  response  and  design  considerations  may  be 
expressed  in  tenns  of  three  basic  parameters:  P  the  proportion  of  time 
in  turbulence,  ac  the  gust  severity,,  and  a  shape  parameter  which  defines 
the  generalized  exceedance  curves.  Over  a  dozen  different  families  cf  I 
theoretical  exceedance  curves  are  generated.  The  basic  response  and  en-l 
vironmental  parameters  that  are  of  concern  in  design  are  discussed,  and 
the  use  of  composite  values  of  these  parameters  as  might  be  involved,  in 
mission  considerations  is  shown.  The  composite  gust  intensity  value  ac, 
and  the  related  scale  value,  L,  still  represent  unsettled  questions.  A 
resexamination  of  some  previous  airline  operational  data  In  terms  of 
generalized  exceedance  curves  is  included.  Recommendations  on  four  spe¬ 
cific  design  procedures  are  given.  One  of  the  design  procedures,  based 
primarily  on  the  response  parameters  A  and  Nq,  incorporates  the  results  | 
of  computational  studies  that  were  performed  on  certain  existing  air-  1 
craft  as  a  means  for  establishing  design  boundaries.  The  problem  of  de-| 
termining  the  probability  of  exceeding  given  load  levels  in  flights  of  | 
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13.  Abstract  (Cont’d) 


specified  duration  is  investigated.  Aread 
and  parameters  which  are  considered  to 
have  weakness  or  uncertainty  are  indicated 
and  recommendations  are  accordingly  made 
for  appropriate  future  research  effort. 
Consideration  of  the  effect  of  filtering 
on  deduced  scale  value  is  included  in  an 
Appendix. 


